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EXECUTIVE  SUMMARY 


This  dociunent  presents  an  annual  report  to  the  Office  of  Naval  Research  for  a  research  program 
entitled  "A  Theoretical  Search  For  Supervelocity  Semiconductors”.  This  program  has  been  fimded  by 
ONR  since  1974  in  the  Department  of  Electrical  and  Computer  Engineering  at  N.C.  State  University. 
The  research  has  resulted  in  more  than  80  refereed  publications  and  numerous  conference  presentations 
from  its  inception.  Major  contributions  to  the  field  of  hot  electron  transport  and  semiconductor  device 
modeling  have  been  achieved,  new  computatitmal  methods  have  been  developed  (e.g.  path  integral 
Monte  Carlo  techniques),  and  the  work  has  helped  stimulate  commercial  ventures  in  the  ai^lications  of 
quaternary  semiconductor  materials  to  electronic  and  optical  devices.  In  addition,  there  have  been 
twenty  eight  Ph.D.  and  M.S.  students  who  have  received  degrees  at  N.C.  State  University  with 
research  support  from  this  contract  Three  visiting  faculty  members  from  Jiq)an  came  to  the  University 
to  work  with  the  faculty  investigators  supported  under  this  ONR  contract  during  the  1979-1983  time 
period.  A  visiting  professor  from  the  French  CNRS  Microstructures  and  Microelectronics  Laboratory 
in  Bagneux  (near  Paris)  spent  a  sabbatical  year  at  N.C.  State  during  1988-89,  and  he  devoted  full-time 
working  on  this  program  at  no  cost  to  ONR.  During  the  current  funding  period,  a  visiting  scholar 
from  China  is  a  member  of  our  research  group  working  on  projects  which  directly  impact  this  ONR 
program.  This  researcher  will  be  appointed  as  a  Visiting  Assistant  Professor  in  August,  1992,  and  will 
spend  approximately  50%  of  his  time  working  on  this  project  His  work  will  be  devoted  to  advanced 
device  simulation  and  physics.  In  particular,  he  will  develop  new  concepts  for  the  real  space  transfer 
logic  transistor  and  other  charge  injection  transistor  structures.  Finally,  we  have  hired  a  new  post¬ 
doctoral  research  associate  from  the  University  of  Michigan  (Prof.  J.  Singh’s  group)  and  he  >viU  join 
us  to  work  on  transport  theory  in  various  quantum-based  devices. 

This  initial  phases  of  this  work  centered  around  the  development  of  Monte  Carlo  simulation  tech¬ 
niques  which  allow  the  study  of  detailed  physics  of  hot  electron  transport  in  a  variety  of  compoimd 


semiconductor  materials.  The  original  emphases  were  concerned  with  electronic  materials  phenomena. 
Later  woiic  considered  the  utilization  of  these  materials  in  realistic  device  structures  where  i^ysical 
boundary  conditions  must  be  imposed  on  the  carrier  transport  More  recently,  the  woik  has  focused  on 
the  domain  of  ultra-small  materials  and  device  i^enomena  where  microscopic  non-local  transient 
effects  such  as  velocity  overshoot  ballistic  and  neaiiy-ballistic  transpoit,  and  quantum  transport 
become  important  or  dominant  During  the  past  five  years  we  have  researched  die  applications  of  the 
Feynman  "integral  over  paths"  iqiproach  to  quantum  transpoit  and  identified  numerical  limitations  to  its 
practical  application.  Also  we  have  emfdiasized  the  study  of  hot  electron  effects  in  new  device  struc¬ 
tures,  such  as  the  hot  electron  spectrometer,  heterojunction  bipolar  transistor,  small  dimension  metal- 
semiconductor-metal  photodetectors,  delta-doped  high  electron  mobility  transistors,  and  real  space 
transfer  logic  transistors.  During  this  time,  we  have  been  exploring  some  new  sqiproaches  to  device 
modeling  which  combine  the  Mrnite  Cario  method  with  the  method  of  moments  of  the  Boltzmann  tran¬ 
sport  equation  (hydrodynamic  transport  model)  for  studying  specific  device  structures,  such  as  small- 
dimension  n^-n-n^  majority  carrier  devices  and  the  high  electron  mobility  transistor.  In  addition,  we 
have  incorporated  quantum  correction  terms  into  the  hydrodynamic  model  and  applied  this  model  to 
resonant  tunneling  structures.  Recently,  we  have  ^>plied  a  new  ansatz  distribution  function  as  a  consti¬ 
tutive  relation  to  close  the  moment  equations  in  the  hydrodynamic  transport  model.  Initial  results  of 
this  approach  have  been  physically  satisfying  and  computationally  promising.  A  comprehensive  formu¬ 
lation  of  this  new  model  has  been  published  and  has  attracted  a  lot  of  attention. 

During  the  next  year,  new  research  in  four  general  areas  will  advance  our  basic  research  direc¬ 
tions  encompassing  the  study  of  hot  electron  transpoit  in  materials  and  devices.  These  four  areas 
included  1)  the  study  of  quantum  transpoit  in  mesoscopic  structures,  nanostructures,  and  related  ultra- 
small  heterobarrier  device  structures  using  several  physical  models  and  computational  approadies,  2) 
new  approaches  to  the  merging  of  Monte  Carlo  methods  with  moment  equatitm  methods  with  resulting 
improvements  to  the  hydrodynamic  transport  model,  3)  Monte  Carlo  simulations  in  order  to  study 


pseudomorphic  devices,  real  space  transfer  structures,  and  submicron  MOSFETs,  and  4)  search  for  new 
quantum-based  device  concepts  for  ultra-fast,  ultra-dense  applications.  The  possibility  of  expanding 
our  experience  to  optical  and  optoelectronic  device  modeling  will  be  seriously  explored  as  well,  as  we 
define  the  role  of  our  new  personnel  which  will  be  added  to  our  group  in  Fall  1992. 
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1.0  INTRODUCTION 

In  October,  1974,  Oie  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research  program 
in  the  Department  of  Electrical  and  Computer  Engineering  at  North  Carolina  State  University.  The 
general  goal  of  this  research  program  has  been  and  continues  to  be  tiie  investigation  of  high-speed  car¬ 
rier  transport  in  m-V  compound  semiconductors,  III-V  alloy  materials,  other  advanced  electronic  and 
optical  materials,  and  novel  device  structures  which  utilize  these  materials.  Four  faculty  members  at 
N.  C.  State  have  been  primarily  involved  and  supported  by  this  project  Currently,  Profs.  M.  A. 
Littlejohn  and  K.  W.  Kim  serve  as  the  co-principal  investigators  on  this  projea.  This  research  pro¬ 
gram  has  made  significant  contributions  to  the  uiKlerstanding  and  knowledge  base  of  hot  electron  tran¬ 
sport  in  materials  and  devices.  It  has  provided  scioitific  guidance  to  the  U.  a.  Navy  in  the  lurmulation 
of  a  part  of  its  basic  and  af^ied  research  program.  Numerical  concepts  develcqied  under  this  project 
have  been  transferred  to  other  Universities,  including  the  University  of  Illinois,  and  we  continue  these 
inter-institutional  collaborations.  During  the  past  year,  we  have  maintained  close  ctxitact  with  the 
research  program  of  Prof.  Hess  at  the  University  of  Illinois.  Also,  collaborations  with  Dr.  H.  L.  Gru- 
bin  of  Scientific  Research  Associates,  Inc.  and  Dr.  G.  J.  lafiate  and  Dr.  M.  A.  Stroscio  of  tiie  U.  S. 
Army  Research  Office  in  the  area  of  quantum  transport  in  semiconductor  devices  have  been 
strengthened.  Dr.  lafrate.  Dr.  Stroscio,  and  Dr.  Giubin  are  Adjunct  Professors  in  the  ECE  Department 
at  North  Carolina  State  Uiuversity  and  we  share  graduate  students  and  post-doctoral  research  associates 
on  joint  projects.  In  addition  to  those  directions,  our  research  results  have  helped  stimulate  commercial 
ventures,  particularly  in  the  develoinnent  of  GalnAsP-based  materials  and  devices.  To-date,  this  pro¬ 
gram  has  resulted  in  79  refereed  publications  in  the  literature,  1  additional  manuscript  is  currently  in 
press,  and  6  manuscriiHs  have  been  submitted  or  are  in  preparation  for  submission  to  the  technical 
literature.  A  listing  of  these  publications  is  given  in  Appendix  A.  In  addition,  numerous  invited  talks 
and  presentations  have  been  given  at  conferences  and  workshops  throughout  the  United  States  and  in 
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Other  countries.  The  program  has  contributed  significantly  to  the  educational  program  at  N.  C.  State 
University  with  more  than  thirty  Ph.D.,  M.S.  and  undergraduate  students  having  received  support  under 
this  ONR  contract.  Currently,  three  Ph.D.  students  are  working  toward  their  degrees  on  this  project 
Two  of  these  Ph.D.  students  are  U.S.  Gtizens.  A  visiting  scholar  from  China,  who  will  be  appointed 
as  a  Visiting  Assistant  Professor  in  August,  1992,  is  also  working  on  research  related  to  this  program. 
We  have  hired  another  post-doctoral  research  associate  who  will  extend  our  cs^rabilities  into  the  area  of 
band  structure  effects,  quantum  transport,  and  many-body  i^ysics. 

This  ONR  program  has  been  efficient  and  productive.  The  high  quality  of  this  research  will  con¬ 
tinue  to  be  maintained  in  the  future.  This  report  summarizes  the  progress  and  accomplishments  made 
during  the  past  contract  period. 
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2.0  RESEARCH  RESULTS 


2.1  Background 

Since  its  inception,  semiconductor  technology  has  been  stimulated  by  requirements  for  electronic 
systems  with  ever-increasing  capabilities  to  process  information  faster,  more  functionally  and  more 
efficiently.  These  requirements  have  motivated  the  scaling  down  of  integrated  circuit  (IQ  device 
dimensions  into  the  sulMnicrai  Gcss  than  ten  thousand  angstroms)  and  ultrasubmicron  Gess  than  one 
thousand  angstroms)  regions.  Today,  we  have  entered  an  era  where  nanostructure  physics  and  fabrica¬ 
tion  motivate  our  research  efforts  in  semiconductor  electronics  [1].  As  fabrication  technology  has 
allowed  such  devices  to  be  realized,  many  new  and  fundamental  questions  have  emerged  concerning 
the  underlying  physics  of  small  (atomic  level)  dimensions  in  semiconductor  devices.  Important  issues 
now  under  consideration  for  ultrasubmicron  devices  include  nonequilibrium  transport  dealing  with  such 
topics  as  quasi-ballistic  transport,  overshoot  i^nomena  and  quantum  transport  A  great  deal  of  pro¬ 
gress  has  been  achieved  in  our  understanding  of  these  important  device  effects,  although  major  work 
remains  to  be  done  as  our  ability  to  fabricate  very  smaU  electronic  device  structures  continues  to 
expand  and  mature  [2,3]. 

The  ability  to  fabricate  small  devices  has  been  continually  refined  over  the  last  ten  years  dirough 
impressive  improvements  in  materials  growth  technologies.  Molecular  beam  epitaxy  (MBE),  metalor- 
ganic  chemical  vapor  deposition  (OMCVD),  and  atomic  layer  epitaxy  (ALE)  have  provided  die  ability 
to  fabricate  a  wide  variety  of  materials  and  heterostructure  combinations  with  near  perfect  interfaces, 
doping  control  and  compositional  uniformity  with  atomic  level  dimensions.  The  development  of  ALE 
may  very  well  prove  to  be  the  ultimate  growth  technology  since  it  allows  the  deposition  of  one  mono- 
layer  of  device  quality  materials  through  a  controllable,  self-limiting  mechanism,  and  is  especially  use¬ 
ful  for  the  deposition  of  heterojunctions  [4].  The  ability  to  grow  layers  with  dimensions  of  a  few 
angstroms  opens  the  domain  of  quantum  transport  to  experimental  study  and  verification.  Thus,  topics 


4 


resulting  from  size  quantization  in  condensed  matter  must  be  investigated  from  theoretical  viewpoints 
with  tools  which  are  either  partially  developed  or  through  the  development  of  new  tools  which  are  not 
now  available.  Quantization  effects  arising  horn  geometrical  size  constraints,  proximity  effects  result¬ 
ing  from  closely  packed  arrays  of  devices,  and  general  solid-state  considerations  not  heretofore  con¬ 
sidered  questionable  (effective  mass  approximation,  the  role  of  contacts  and  the  like)  must  be 
addressed  from  a  fundamoital  point  of  view.  Moreover,  from  a  device  physics  point-of-view.  it  is 
desirable  to  have  a  microscopic  description  of  the  physics  of  smaU  dimensions  which  is  amenable  to 
phenomenological  treatment,  so  that  its  properties  can  be  meaningfully  incoiporated  into  futuristic  dev¬ 
ice  concepts  and  simulations. 

Theoretical  methods  to  address  carrier  transport  have  also  progressed  rapidly  over  the  last  ten- 
fifteen  years,  in  a  manner  similar  to  research  in  semiconductor  thin  film  epitaxial  growth  technology. 
In  fact,  this  is  a  natural  progression  in  many  ways  and  is  to  be  expected.  The  progress  achieved  in 
materials  growth  of  structures  with  quantum  dimensions  dictates  that  new  jq>proaches  be  developed  and 
refined  to  study  quantum  transport  phenomena  and  the  physics  of  small  dimensions.  However,  a 
significant  change  in  direction  is  now  warranted.  Past  transport  theory  and  device  modeling 
approaches  have  relied  on  particle  or  quasi-particle  approaches  where  the  electrons  are  treated  as  rigid 
mobile  entities  which  undergo  interactions  with  the  transport  medium.  The  treatment  of  the  interaction 
often  involves  wave  concepts.  However,  the  model  is  basically  a  particle  model.  In  the  current  regime 
of  quantum  transport,  we  may  no  longer  be  able  to  consider  the  carriers  as  particles.  It  is  quite  likely 
that  their  physical  behavior  will  be  governed  either  partially  or  completely  by  wave  phenomena. 

Quasi-particle  methods  are  an  attempt  to  model  structures  with  quantum  dimensions  which  retain 
as  much  of  the  classical  formalism  as  possible  in  order  to  be  able  to  express  results  in  terms  of  param¬ 
eters  which  are  of  the  greatest  experimental  interest,  such  as  carrier  velocity  and  diftusion  constant 
This  is  a  flexible  approach.  However,  much  care  is  required  to  ensure  that  aU  important  effects  are 


5 


properly  included  because  of  the  approximations  involved  in  the  formulations.  On  the  other  hand, 
more  fundamental  quantum  approaches,  such  as  operator-eigenfunction  methods,  adhere  closely  to  the 
actual  quantum  states  present  in  the  device  stnicture  when  scattering  is  not  included.  Scattering 
processes  (dissipation)  can  be  added  by  using  perturbation  theory  from  quantum  mechanics.  These 
techniques  can  obtain  the  greatest  sensitivity  to  the  resulting  carrier  confinement  and  die  lattice  poten¬ 
tials.  However,  they  are  relatively  inflexible  in  studying  non-linear  dynamical  properties  in  the  pres¬ 
ence  of  strong  dissipation,  such  as  is  present  in  tte  electron-phonon  interaction  at  high  electric  fields. 

Another  approach  to  quantum  transport  relies  on  the  "integral  over  paths"  method,  originally  pro¬ 
posed  by  Dirac  and  formulated  by  Feynman  [5).  Practical  path  integral  methods  for  the  study  of  small 
devices  rely  on  an  influence  functional  technique  in  which  the  source  of  the  dissipation  has  been 
integrated  over  all  phonon  modes.  This  results  in  a  model  influence  functional  where  the  phonon¬ 
scattering  dissipation  can  be  represented  as  an  interaction  with  a  collection  of  harmonic  oscillator 
modes  in  which  the  translational  invariance  of  the  carriers  is  preserved.  The  resulting  model  includes 
constant  or  oscillatory  electric  and  magnetic  fields,  carrier  screening,  scattering  and  dissipation,  carrier 
confinement,  background  temperature  and  initial  conditions  can  be  dealt  with  as  readily  as  for  a  free 
particle  [6].  The  path  integral  method  can  be  compared  and  contrasted  to  other  methods  now  being 
studied  and  supported  by  ONR,  such  as  the  ai^lication  of  the  density  matrix  formalism  and  Wigner 
distribution  function  s^proaches  including  the  use  of  moment  equations  [7].  Based  on  our  progress  in 
developing  the  path  integral  approach  during  the  past  five  years,  we  believe  that  path  integral  methods 
will  play  an  important  role  in  developing  an  understanding  of  quantum  transport  in  reduced  geometries 
(such  as  mesoscopic  devices  [8]).  However,  because  of  its  intensive  computation  requirements,  the 
path  integral  method  is  very  limited  in  its  application  to  practical  devices.  CTurrently,  we  are  aware  of 
only  one  other  U.S.  university  program  directed  toward  iq)plying  the  path  integral  method  to  semicon¬ 
ductor  devices.  This  is  the  program  at  the  University  of  Illinois  under  the  direction  of  Prof.  Karl  Hess. 
Our  efforts  in  the  development  of  this  method  are  complimentary  to  the  one  at  Illinois,  and  we  have 
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collaborated  with  the  Illinois  personnel  through  the  use  of  the  NSF  Superc(»nputing  Center  and  the 
National  Center  for  Computational  Electronics  located  there.  One  of  our  former  students,  whose  Ph.D. 
dissertation  on  path  integral  methods  was  supported  by  this  ONR  contract,  joined  Prof  Hess’  group  in 
January,  1990  as  a  Visiting  Research  Assistant  Professor.  We  plan  to  maintain  contact  with  Prof. 
Hess’  group  as  we  evaluate  our  path  integral  research  directions.  We  will  continue  this  work  if  we  can 
hire  an  appropriate  Ph.D.  candidate. 

The  work  at  N.  C.  State  University  supported  by  the  ONR  under  die  current  contract  has  pro¬ 
gressed  over  the  past  fifteen  years  from  the  realm  of  particle  models  to  quasi-particle  models  to  quan¬ 
tum  transport  models.  We  have  relied  primarily  on  the  Monte  Carlo  method  to  study  and  solve  tran¬ 
sport  problems  in  III-V  compound  semiconductors.  Our  techniques  for  modeling  materials  physics  and 
device  phenomena  extend  to  device  dimensions  around  one  thousand  angstroms  [10,11],  and  we  are 
confident  in  these  models  for  predicting  steady  state  and  transient  device  effects  down  to  these  dimen¬ 
sions.  However,  in  order  to  remain  in  the  forefront  of  trartsport  j^ysics  and  device  research,  we  must 
continue  and  increase  our  progress  into  the  realm  of  dimensions  where  physical  effects  can  be  studied 
in  device  regions  with  dimensions  less  than  one  thousand  angstroms.  We  will  continue  to  explore  the 
Monte  Carlo  ttanspcrr.  method  for  study  of  novel  device  structures  where  theoretical  underpinning  is 
required.  One  novel  class  of  structures  which  will  be  studied  in  detail  by  the  Monte  Carlo  method  is 
based  on  the  rapidly-emerging  pseudomorphic  or  strained-layer  devices.  These  structures  permit 
extended  compositional  ranges  and,  thus,  have  a  number  of  potential  advantages  such  as  higher  tran¬ 
sconductance  and  channel  carrier  density.  Recent  studies  on  the  strain-induced  piezoelectric  fields  have 
op)ened  yet  another  possibility  in  realizing  ultra-fast  switching  devices,  resulting  in  an  increased  impor¬ 
tance  on  the  study  of  carrier  transport  in  these  devices.  Other  novel  devices  which  will  be  analyzed 
include  delta-doped  field-effect  transistors  and  real-space  transfer  devices,  based  on  studies  which  have 
begun  during  the  past  two  years.  In  addition,  we  will  incorporate  Monte  Carlo  methods  into  the 
quasi-particle  approach  based  on  moments  of  the  Boltzmann  transport  equation  with  quantum 
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mechanical  corrections.  Preliminary  results  show  the  utility  of  this  approach  for  modeling  microwave 
and  millimeter  wave  devices  which  are  important  to  the  DoD’s  MMIC  program. 

Finally,  we  want  to  re-iterate  the  important  interactions  which  have  developed  between  our 
research  group  and  other  researchers  during  the  last  year.  These  inten.vtions  increase  the  impact  of 
this  program  on  the  field  of  semiconductor  device  physics  and  allow  our  research  efforts  to  be  far  more 
productive  through  increased  intellectual  efforts  and  enhanced  facilities  and  resources.  During  the 
1991-1992  funding  period,  we  have  maintained  collaborations  with  Prof.  K.  Hess  at  the  University  of 
Illinois,  Dr.  H.  Grubin  of  Scientific  Research  Associates,  and  Drs.  M.  Stroscio  and  G.  lafrate  of  the 
U.  S.  Army  Research  Office.  These  research  collaborations  have  resulted  in  several  joint  publications 
based  on  mutual  research  interests,  expertise  and  c^abilities.  There  have  been  several  visits  between 
these  laboratories  and  the  logistics  for  increased  collaborations  are  excellent. 

Thus,  we  believe  that  we  have  a  well-rounded  approach  for  this  ONR  program  which  will  lead  to 
new  innovations,  and  practical  developments  in  semiconductor  device  physics,  simulation  and  model¬ 
ing. 

2.2  Summary  of  Research  Results 

This  section  provides  the  current  status  of  our  research  efforts  in  specific  topical  areas  we  have 
worked  during  the  past  year. 

2.2.1  Hydrodynamic  hot-electron  transport  model 

The  primary  goal  of  this  research  is  to  develop  and  implement  an  advanced  hot-electron  transport 
model  suitable  for  accurate  and  efficient  study  of  electron  dynamics  in  both  submicron  and  ultrasmall 
semiconductor  device  structures.  The  model  is  required  to  accurately  reflect  hot-electron  effects  such 
as  velocity  overshoot  and  nonstationary  effects  due  to  spatial  nonuniformities  in  average  electron  con¬ 
centration,  average  electron  velocity  and  average  electron  energy.  Also,  the  model  is  required  to  accu- 


8 


rately  describe  the  streaming  motion  and  dissipation  of  electrons  in  multiple  noiqjarabolic  conduction 
bands  (as  are  present  in  especially  all  III-V  compound  semiconductors).  In  addition,  an  electron  tran¬ 
sport  model  is  desired  which  can  be  extended  to  the  study  of  quantum  systems  (i.c.,  feamre  sizes  less 
than  1000  A)  where  wave-like  characteristics  (e.g.,  tunneling  and  interference  effects)  are  important. 

Over  the  past  decade,  much  effort  has  been  focused  on  the  development  of  semi-classical  electron 
transport  models  which  accurately  reflect  hot-electron  effects  in  submicron  semiconducting  device 
structures  [11].  The  incentive  is  that  nonequilibrium  and/or  nonstationary  effects  may  be  employed  to 
achieve  superior  electron  device  performarKe.  One  important  example  of  a  potentially  useful  hot- 
electron  effect  is  velocity  overshoot  [12].  In  velocity  overshoot,  the  difference  in  momentum  and 
energy  relaxation  times  can  result  in  an  initial  temporal  drift  velocity,  under  constant  electric  field  con¬ 
ditions,  which  overshoots  the  corresponding  time  steady-state  (stationary)  value.  This  dynamical 
response  in  electron  motion  is  due  to  the  difference  between  the  displacement  of  the  election  distribu¬ 
tion  in  momentum  space  ajtd  the  displacement  in  energy  space  (i.e.,  the  form  of  the  nonequilibrium 
distribution  function).  For  this  type  of  election  motion,  one  cannot  use  simple  models  based  upon  sta¬ 
tionary  relationships  (which  ignore  important  features  of  the  distribution)  to  describe  the  transport  phy¬ 
sics.  Furthermore,  deficiencies  in  these  near-equilibrium  models  may  also  fail  to  accurately  describe 
electron  transport  in  submicron  device  structures  even  in  time  steady-state.  This  can  occur  when  the 
election  gas  is  far  from  equilibrium  and  large  spatial  gradients  :n  the  electron  gas  parameters  (electron 
density,  average  electroi.  /elocity,  average  election  energy,  etc.)  exist. 

Since  classical  (drift-diffusion)  models  assume  local  stationary  relationships  (average  velocity 
depends  on  the  local  electric  field)  and  are  based  upon  near  equilibrium  derivations  [12],  physically 
superior  models  must  be  used  to  fiiUy  and  effectively  study  important  hot-electron  effects.  At  present, 
the  problem  of  treating  transport  in  submicron  structures  is  typically  attacked  using  one  of  two  distinct 
methods.  The  first  and  most  accurate  approach  has  been  to  apply  Monte  Carlo  simulation  methods 
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[13].  The  Monte  Carlo  method  is  a  powerful  statistics-based  numerical  technique  for  investigating 
electron  transport  physics  in  semiconductors  and  semiconductor  devices  [13].  While  the  Monte  Carlo 
method  provides  a  relatively  simple  and  accurate  indirect  approach  to  determine  particle  distribution 
functions  and  transport  parameters  resulting  from  the  Boltzmarm  transport  equation  subject  to  complex 
scattering  mechanisms,  the  approach  does  have  important  disadvantages.  Specifically,  the  technique 
often  requires  extensive  computation  time  to  arrive  at  adequate  statistics  and  does  not  naturally  possess 
macroscopic  terms  for  physical  interpretation. 

The  alternative  (second)  approach  to  study  electron  device  physics  more  accurately  and  efficiently 
is  to  develop  macroscopic  models  which  predict  average  transport  quantities  as  opposed  to  the  full 
electron  distribution  function.  In  this  approach,  the  description  of  electron  dynamics  is  reduced  to 
hydrodynamic-like  conservation  equations  for  electron  density,  average  electron  momentum,  average 
electron  energy,  etc  [14].  The  hydrodynamic  models  are  derived  from  moments  of  the  Boltzmann  tran¬ 
sport  equation  which  usually  include  many  simplifying  assumptions  (i.e.,  usually  neglect  conduction 
band  nonparabolicity,  treat  multiple  conduction  band  valley  effects  indirectly  and  directly  or  indirectly 
make  incorrect  assumptions  about  the  form  of  the  electron  distribution  function)  about  the  electron 
transport  physics.  However,  if  advanced  hydrodynamic  models,  which  include  hot  electron  eflects 
such  as  multiple  nonparabolic  conduction  bands,  are  developed  and  combined  with  Monte  Carlo  tech¬ 
niques  this  is  an  excellent  overall  approach.  A  limited  number  of  Monte  Cario  simulations  can  then  be 
utilized  to  verify  the  accuracy  of  the  simplified  models  and  to  supply  approximate  parametric  dissipa¬ 
tion  data  for  the  particular  materials  and/or  devices  under  consideration.  These  types  of  advanced 
hydrodynamic  models,  together  with  the  Monte  Cario  method,  represent  a  simplified  mathematical 
approach  for  accurate  and  efficient  study  of  realistic  transport  physics. 

Finally,  very  few  hydrodynamic  (macroscoinc  or  average  ensemble-electron)  models  have  been 
implemented  to  accurately  represent  quantum  effects  (i.e.,  tunneling  and  interference)  present  in 
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ultrasmall  electron  structures.  One  method  to  accurately  model  structures  with  spatial  dimensions  on 
the  scale  of  the  de  Broglie  wavelength  (i.e.,  in  the  quantum  regime)  is  to  utilize  quantum  hydrodynam¬ 
ics  [15,16].  In  this  method,  a  set  of  quantum  moment  equations  are  derived  by  taking  moments  [17] 
of  the  Wigner-Boltzmann  equation  and  closure  is  obtained  using  an  tq>proximate  nonequilibrium  form 
of  the  zero-current  (equilibnum)  Wigner  distribution  function.  This  procedure  results  in  a  set  of 
quantum-corrected  hydrodynamic  transport  equations  which  have  been  demonstrated  to  predict  quantum 
pheimmena  [15,16]. 

During  the  past  year,  advanced  forms  of  the  hydrodynamic  transport  model  (HTM)  were 
developed  and  utilized  to  study  both  submicron  and  ultrasmall  semiconducting  device  structures.  First 
an  advanced  semi-classical  hydrodynamic  transport  model  was  developed  and  implemented  to  accu¬ 
rately  predict  electron  dynamics  in  multiple  nonparabolic  energy  bands.  In  addition,  a  new  quantum 
mechanical  formulation  of  the  first  three  moment  equations  (i.e.,  quantum  hydrodynamic  equations  for 
electron  density,  average  velocity  and  average  energy)  was  investigated  and  implemented  to  predict 
non-classical  timneling  in  a  double-barrier  heterostructure.  The  development  of  the  advanced  model(s) 
was  achieved  through  several  separate  theoretical  investigations. 

First,  to  determine  a  method  to  include  the  effects  of  conduction  band  nonparabolicity  into  the 
macroscopic  streaming  (collision  independent)  of  the  semi-classical  HTM,  the  fundamental  physics  of 
electron  transport  were  investigated.  An  alternate  formalism  was  developed  for  deriving  a  new  nonpar¬ 
abolic  hydrodynamic  transport  model  adequate  for  modeling  hot  electron  transport  in  submicron  sem¬ 
iconductor  devices.  Specifically,  the  model  equations  were  developed  by  applying  a  unique  set  of 
moment  operators  (Oo(k)=l;  a  constant,  0|(k>=u(k);  the  nonparabolic  velocity,  4>2(k)= 
[m(k)/2]u(k)  u(k>=E(k);  an  ^roximation  to  the  particle  energy)  to  the  collisionless  Boltzmann  tran¬ 
sport  equation.  The  operators,  used  in  this  analysis,  represent  a  significant  improvement  in  the 
development  of  the  hydrodynamic  transport  model  since  they  lead  to  a  more  compact  mathematical 
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form.  This  resulting  form  can  be  manipulated  more  easily  and  reveals  clearly  the  allowable  simplify¬ 
ing  approximations.  An  intuition-based  distribution  function  was  proposed  as  a  constitutive  (essential) 
relation  to  close  the  moment  equations.  This  ansatz  nonparabolic  distribution  function,  which  is 
derived  in  velocity  space  for  mathematical  convenience  and  approximates  some  of  the  characteristics  of 
the  tnre  nonparabolic  distribution  function,  was  then  analyzed  and  used  to  develop  a  set  of  transport 
parameters.  After  devdopment,  the  results  from  two  independent  Monte  Carlo  transport  models  were 
used  to  evalurre  die  nonparabolic  transport  parameters.  These  models  were  used  to  investigate  electron 
transport  in  both  stationary  (uniform  electric  field)  and  nonstationary  (ballistic  diode)  envirorunents. 
The  results  of  these  studies  show  that  there  is  excellent  agreement  between  this  new  nonparabolic  tran¬ 
sport  model  and  the  Monte  Carlo  calculations.  Hence,  this  new  nonparabolic  transport  model  presents 
a  viable  alternative  to  studying  realistic  electronic  device  structures  operating  under  high  bias  condi¬ 
tions. 

In  the  previous  analysis,  specific  macroscopic  collision  terms  are  omitted.  This  is  done  to  permit 
a  concise  focus  on  the  streaming  (collision  independent)  terms  that  evolve  from  the  application  of  non¬ 
parabolic  conduction  bands.  The  collision  terms,  which  are  almost  always  treated  phenomenologically, 
has  been  developed  separately  for  a  nonparabolic  multi-valley  system  (GaAs).  SpecificaUy,  multi¬ 
valley  dissipation  is  included  in  the  model  using  an  ensemble  relaxation  time  approximation.  The 
numerous  intravalley  and  intervalley  relaxation  times  were  calculated  using  the  Monte  Carlo  Method. 

In  order  to  extend  the  modeling  approach  to  the  quantum  device  regime,  a  technique  was  investi¬ 
gated  to  develop  quantum  corrections  for  the  hydrodynamic  equations.  The  set  of  quantum-corrected 
transport  equations  which  we  have  developed  arxl  investigated  were  based  upon  quantum  moment 
equations  of  Grubin  and  Kteskovsky  [15].  Grubin  and  Kreskovsky  previously  constructed  a  set  of 
quantum  balance  (hydrodynamic)  equations  applicable  for  electron  transport  in  mesoscopic  structures. 
Their  model  was  derived  from  the  general  moment  equations  of  Stroscio  [17]  using  the  displaced  none- 
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quilibrium  Wigner  distribution  function  of  Ancona  and  lafrate  [18].  These  resulting  quantum  hydro- 
dynamic  equations  incorporate  quantum  effects  tiuough  the  introduction  of  terms  (corrections)  with 
density  gradient  dependence. 

Both  the  advanced  semi-classical  and  the  quantum-corrected  hydrodynamic  transport  models 
(HTMs)  were  used  to  study  electron  transport  in  novel  semicrxiductor  device  structures.  A  general 
numerical  algorithm  was  implemoited  to  solve  the  resulting  system  of  nonlinear  differential  equations 
[from  both  physical  model(s)]  for  arbitrary  one-space-dimensional  n-type  structures.  The  numerical 
algorithm  utilizes  a  local  simultaneous-coupled  (Newton)  iteration  approach  in  tandem  with  a  global 
physics-based  continuation.  Simulation  results  were  generated  for  a  single-valley  (T)  HTM,  with  quan¬ 
tum  corrections,  ^lied  to  a  double-barrier  heterostnicture  with  Al^Caj.^As  barriers.  This  study  was 
the  first  to  use  the  three  quantum  hydrodynamic  (QHD)  equations  to  study  transport  in  heterostructure 
devices.  The  key  additional  feature  of  the  QHD  equations  is  the  incorporation  of  the  effects  of  density 
gradients  through  the  quantum  potential.  Insofar  as  a  classical  solution  does  rx>t  exist  in  the  presence 
of  barriers,  Q  (i.e.,  one  of  the  major  quantum  correction  terms)  must  at  least  partially  cancel  the  effects 
of  the  barriers  and  permit  carrier  transport.  Indeed,  as  is  shown  in  Fig.  1,  Q/3  does  approximately  bal¬ 
ance  the  barrier  potential  which  is  a  result  consistent  with  Ref.  15.  In  this  woric,  the  quantum 
corrected  HTM  was  used  to  demonstrate  that  the  hydrodynamic  approach  can  be  used  to  study 
ultrasmall  electron  devices. 

In  another  simulation  study,  a  simplified  form  of  the  multi-valley  hydrodynamic  transport  model 
was  used  to  investigate  the  feasibility  of  utilizing  the  velocity  overshoot  effect  in  reducing  electron 
transit  time  through  submicron  GaAs  electronic  device  structures.  In  this  initial  investigation  with  the 
new  multi-valley  model,  which  incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduc¬ 
tion  band  (T,  L,  and  X)  and  allows  for  elearon  exchange  through  energy  dependent  relaxation  parame¬ 
ters,  the  nonparabolic  correction  terms  were  be  suppressed.  This  was  done  to  simplify  the  analysis  aiKi 
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Figure  1:  Conduction  band  (Eg)  and  quantum  potential  (Q)  at  VmppiM  =  8mV.  The  source 
and  collector  regions  are  excluded  from  the  plot. 
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to  provide  a  base  for  future  nonparabolic  studies.  This  new  multi-valley  model  was  afqslied  to  a  novel 
submicron  GaAs  ballistic  diode  stnicture  with  bandgap  engineered  electric  field  spikes  (and/or  doping 
spikes)  in  order  to  reduce  upper  valley  occupancy  and  enhance  device  transit  time  using  the  velocity 
overshoot  effect  The  general  results  of  the  study  (summarized  in  Fig.  2)  indicate  that,  while  velocity 
speed-up  is  miiumal  in  die  neigtdwrhood  of  die  active  (low-doped  or  high-electric  field)  region,  velo¬ 
city  overshoot  can  be  extended  into  the  collector  Oow-electric  field)  region  with  significant  transit  time 
enhancement  Since  these  results  are  consistent  with  other  independent  investigations  using  more  accu¬ 
rate  (Monte  Cario)  techniques,  this  work  has  proven  that  HTMs  are  a  viable  alternative  to  these  more 
computationally  intensive  methods. 

We  have  also  determined  the  impact  of  k-space  transfer  (i.e.,  intervalley  transfer)  and  band  non- 
parabolicity  on  electron  were  determined.  In  this  work,  variations  of  the  multi-valley  (T,  L,  and  X) 
nonparabolic  hydrodynamic  transport  model  was  used  to  study  a  submicron  GaAs  ballistic  diode. 
Numerical  simulations  indicate  that  accurately  including  the  effects  of  nonparabolicity  in  the  streaming 
terms  and  k-space  transfer  in  the  velocity  and  energy  equations  is  very  important  in  correctly  determin¬ 
ing  the  conductance  of  the  device.  The  existence  and  amount  of  negative  difTerential  conductance  was 
determined  to  be  strongly  influenced  by  both  of  these  physical  factors.  Furthermore,  the  sensitivity  of 
device  conductance  to  changes  in  the  thermal  conductivity  is  diminished  significantly  when  nonparabol¬ 
icity  is  accurately  incorporated.  Results  for  the  total  current  density  versus  applied  bias,  for  different 
variations  in  the  model  are  given  in  Fig.  3.  These  simulations  have  demonstrated  the  pronounced 
effects  of  conduction  band  nonparabolicity  and  k-space  transfer  on  electron  transport  in  a  GaAs  ballis¬ 
tic  diode.  Therefore,  the  study  of  nonequilibrium  and  nonstationary  electrcm  i^ysics  should  be  per¬ 
formed  with  this  type  of  advanced  macroscopic  transport  model. 

More  recently,  work  has  been  performed  to  develop  an  extended  Scharfetter-Gummel  current- 
density  discretization  formula  sufficient  for  solving  an  advanced  hydrodynamic  electron  transport  model 
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Figure  2:  Outline  of  traniport  study  and  transit  time  results 
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Figure  3:  Total  electron  current  density  versus  applied  electric  potential  for  a  submicron 
ballistic  diode  (BD  :see  inset).  Bruits  are  for  variations  in  the  hydrodynamic  transport 
model  (HTM  ;  dashed  lines)  and  a  Monte  Carlo  study  (MC  ;  solid  line).  Four  different 
variations  (cases)  in  the  HTM  were  consider.  The  results  for  these  different  cases  are 
designated  by  duhed  lines  with  dash  lengths  that  increase  in  the  following  order.  Case  A 
:  HTM  without  nonparabolicity  (in  the  streaming  terms)  or  full  intervaUey  effects  (in  the 
velocity  or  energy  equation).  Case  B  :  HTM  with  intervalley  effects  included  only.  Case  C  : 
HTM  with  nonparabolicity  effects  included  only.  Case  D  :  HTM  with  nonparabolicity  and 
intervalley  effects. 
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(HTM).  While  this  extended  formulation  properly  incorporates  many  advanced  effects  associated  with 
hot-electron  transport,  its  primary  feature  is  the  accurate  treatment  of  the  convection  term  (or  displace¬ 
ment  energy  term)  which  is  present  in  both  conventional  (i.e.,  drifted  Maxwellian  model)  and  the 
advanced  nonparabolic  model  previously  discussed  in  this  report.  In  general,  since  the  success  of  any 
numerical  simulation  (for  a  wide  variety  of  device  structures  and  operating  conditions)  is  critically 
dependent  on  convergence  rate  and  numerical  stability,  the  choice  of  an  efficient  discretization  scheme 
is  very  important.  Currently,  most  robust  discretization  formulations  for  HTMs  utilize  extensions  of 
the  Scharfetter-Gummel  (SG)  scheme  which  was  found  to  be  stable  in  drift-diffusion  (DD)  model 
simulations.  However,  most  extensions  of  the  SG  formula  for  HTMs  either  ignore  convention  terms  or 
else  assume  a  priori  that  they  are  defined  in  a  previous  step  of  an  iteration  procedure.  A  simple 
analytical  approach  for  directly  including  convection  effects  into  an  extended  SG  formulation  has  been 
established  and  several  versions  of  this  discretization  are  under  investigation.  Presently,  the  stability 
and  efficiency  of  this  new  discretization  formulation  is  still  under  investigation. 

Finally,  work  has  been  performed  to  interrogate  the  influence  of  different  device  structures  and 
boundary  conditions  on  the  resultant  heat  flow  vector  (a  transport  parameter  in  the  HTM)  using  the 
Monte  Carlo  Approach.  Thus  far,  we  have  observed  that  the  sign  of  the  heat  flow  vector  is  a  ftmction 
of  position  in  ballistic  diode  structures  and  that  bourxiary  conditions  do  influence  the  effective  conver¬ 
gence  rate  (i.e.,  simulation  time  to  achieve  correct  symmetry  of  the  electron  distribution)  of  the  Monte 
Carlo  method. 

2.2.2  Modeling  arxl  characterization  of  novel  high  speed  devices 

During  the  past  contract  period,  we  have  studied  electron  transport  properties  and  electrical 
characteristics  for  a  number  of  high  speed  device  stnictures,  including  delta-doped  GaAs  MESFETs, 
AllnAsA3aInAsAnP  HEMTs  with  different  additional  channel  doping  characteristics,  novel  real-space 
transfer  logic  transistor  structures,  and  devices  with  variable  doping  profiles.  The  motivations  of  our 
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research  on  these  stnictures  are  to  provide  in-depth  understanding  of  elecmm  transpoit  in  devices,  to 
investigate  the  possibility  of  realizing  novel  hi^  speed  functional  device  structures,  and  to  provide 
general  design  guidelines  for  optimal  device  operation. 

As  a  promising  high  speed  device,  the  delta-doped  GaAs  MESFET  features  high  current  drive 
capability,  high  transconductance,  improved  threshold  voltage  control,  and  improved  breakdown 
characteristics.  The  study  of  properties  of  delta-doping  in  GaAs  predicts  that  device  performance  of 
delta-doped  GaAs  MESFETs  is  comparable  or  even  better  than  those  of  AlGaAs/GaAs  HEMTs  [19]. 
However,  experimental  results  show  a  wide  diversity  and  fall  short  of  theoretical  predictions.  The 
design  and  optimization  of  device  structure  (especially  delta-doping  profile)  play  an  important  role  in 
guiding  the  improvement  of  device  performance.  Another  high  speed  device,  the  AllnAsAjalnAs/lnP- 
based  HEMT,  has  been  the  focus  of  extensive  experimental  study  due  to  superior  electron  transport 
properties  in  this  lattice-matched  material  system.  While  promising  dc  arxl  microwave  performances 
have  also  been  achieved,  high  output  conductance  has  been  repeatedly  observed  for  devices  with 
channel-length  in  the  deep-submicron  (less  than  0.3  pm)  regime  [20,21].  Physics  insight  leading  to  the 
high  output  conductance  and  potential  solutions  merit  careful  investigation  if  potential  advantages  of 
AlInAsAjalnAs  HEMTs  are  to  be  fully  exploited.  The  third  device  structure  which  we  have  investi¬ 
gated  this  year  belongs  to  a  new  class  of  devices  which  employs  the  concept  of  electron  real-space 
transfer  [22]  in  heterojunction  structures.  The  charge  injection  transistor  (CTHINT)  [23]  and  the 
recently-developed  CHINT  logic  element  [24]  are  representative  devices.  Theoretical  predictions  and 
experimental  evidence  from  CHINT  operation  indicate  the  possibility  for  developing  novel  and 
improved  multi-terminal  real-space  transfer  devices  which  can  realize  multiple  logic  functions.  Our 
last  work  in  the  area  of  device  modeling  deals  with  the  effects  of  electron  transport  in  devices  with 
variable  doping  distributions.  Such  device  stnictures  can  be  realized  by  utilizing  advanced  fabrication 
techniques,  such  as  focused  ion  beams  or  epitaxy  (MBE,  MOCVD,  ALE)  methods.  The  potential 
advantages  of  variable  doping  in  device  channel  or  active  transpon  region  can  best  be  realized  when 
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detailed  electron  transport  properties  in  practical  device  structures  are  understood  and  doping  scheme 
optimization  is  achieved.  In  the  following,  we  outline  the  main  results  of  our  theoretical  study.  More 
detailed  descriptions  for  each  subject  can  be  found  in  the  related  papers  and  presentations. 

A.  Delta-doped  GoAs  MESFET 

Figure  4  gives  a  schematic  illustration  of  the  device  structure  which  we  have  simulated,  along 
with  the  delta-doping  profile.  In  this  work,  we  use  a  two-dimensional,  d^ft-diffusion  model  to  exam¬ 
ine  the  influence  of  key  design  parameters  on  device  operation.  These  parameters  ate  (Fig.  4)  peak 
doping  density,  half  width  at  half  maximum  of  the  doping  profile,  and  distance  from  the  location  of 
peak-doping  to  the  device  surface,  top  layer  background  doping,  and  lateral  feature  sizes.  The  general 
design  considerations  and  trade-offs  are  studied.  In  particular,  our  simulation  results  reveal  that:  1)  A 
steeper  doping  profile  and  a  short  gate-to-channel  distance  are  preferred  for  improvement  in  device  per¬ 
formance.  2)  The  peak  delta-doping  concentration  and  half-width  at  half  maximum  of  delta-doping 
pr  tie  are  found  to  have  a  more  profound  effect  than  does  delta-doping  depth  (dp)  in  enhanced  overall 
device  performance.  And  3)  careful  trade-off  should  be  made  in  the  design  of  top-layer  background 
doping  in  terms  of  current  drive,  transconductance,  and  threshold  voltage  sensitivity  control.  Our 
results  show  that  a  relatively  high  top  layer  background  doping  is  preferred  to  increase  drive  current 
and  reduce  relative  threshold  voltage  shift  (due  to  processing  parameter  fluctuations)  at  a  cost  of  a 
modest  reduction  in  transconductance.  Simulation  results  are  in  close  agreement  with  measurements 
from  experimental  devices  f.?  icated  at  Texas  Instruments  using  atomic  layer  epitaxy  material  prepared 
in  NCSU. 

B.  AlInAsiGalnAs  HEMT 

Figure  5  shows  the  simulated  quarter-micron-gate  AIInAsA3aInAs  HEMT  structure  which 
employs  different  additional  channel  doping  (n^,  i,  and  p)  configurations.  The  motivation  is  to  study 


(a) 


Figure  4:  (a)  Schematic  illustration  of  a  delta-doped  GaAs  MESFET.  The  mesh  distribution  is  illus¬ 
trated  in  the  inset  b)  Delta-doping  profile.  The  solid  line  shows  the  measurement  data  and  the  dashed 
line  shows  die  Gaussian  distribution  used  in  the  simulations. 
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Figure  5:  Schematic  drawing  of  the  layer  structures  for  AlInAs/GalnAs  HEMTs  used  in  this  study. 
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the  effect  of  additional  channel  doping  on  device  transconductance,  output  conductance,  and  high  fre¬ 
quency  performance.  The  AlInAs/GalnAs/InP  material  system  appears  to  be  the  most  promising  candi¬ 
date  for  high  frequency  applications.  However,  for  this  material  system,  high  output  conductance  has 
been  an  major  problem  that  degrades  high  frequency  performance  when  the  HEMT  device  feature  size 
is  reduced  to  the  ultrasubmicron  (less  than  0.3  pm)  regime.  This  has  been  reported  by  several  research 
groups.  Attempts  have  been  made  experimentally  to  reduce  output  conductance,  however,  physical 
insight  remains  an  unknown.  In  this  woik.  we  investigated  the  possible  mechanisms  which  lead  to 
enhanced  or  degraded  device  performance  in  different  char  et-doped  structures  using  two-dimensional, 
self-consistent  ensemble  Monte  Carlo  simulations.  Our  results  indicate  that  conduction  in  the  AlInAs 
buffer  layer  and  reduced  electric  field  shielding  in  n-channel  doped  HEMT  are  the  possible  reasons  for 
increased  output  conductance.  B;  use  of  a  properly  designed  p-doped  channel  layer,  we  observe 
significantly  reduced  output  conductance  and  increased  tra^^^'  .  Juctance  compared  with  those  of  n- 
charmel  '  oped  HEM '  i.  Detailed  analysis  ir^dicates  that  improved  electron  confinement  in  the  device 
channel  anu  increase,  leld  shieding  fcr  HEMT  with  p-channel  doping  improve  overall  device  perfor- 
rr  nee.  Our  Simula  is  are  in  reasonable  agreement  with  measurement  results  from  experimentally 
fabr  .ated  devices  (with  different  additional  channel  doping  (UCSB  and  Hughes  Research  Lab)]. 

C.  Real-space  transfer  logic  transistor 

The  basic  real-space  transfer  logic  transistor  (RSTLT)  is  a  four-terminal  device,  which  has  two 
separate  collectors  (Cl  and  C2)  as  output  terminals  when  connected  through  load  resistors  (R^i  and 
Rc2)  to  the  collector  supply  voltage  (Vc);  a  drain  input  terminal  (D);  and  a  common  source  terminal 
(S)  which  is  at  ground  potential  (Fig.  6).  The  operational  principle  of  the  proposed  RSTLT  is  based 
on  the  localized  nature  of  real-space  transfer  of  hot  electrons  in  such  structures.  Channel  electrons 
gain  energy  from  the  lateral  electric  field  and  real-space  uansfer  occurs  when  the  energy  of  electrons  is 
comparable  to  or  higher  than  the  heterointerface  barrier  height  The  transferred  electrons  will  be  col- 


Figure  6:  Proposed  real-space  transfer  logic  transistor  (RSTLT)  structure. 
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lected  by  either  Cl  or  C2  through  vertical  electric  held.  For  a  given  heterointerface  barrier  height,  the 
spatial  location  of  electron  RST  is  determined  primarily  by  the  bias  conditions.  By  properly  designing 
the  device  configuration  and  selecting  the  operating  voltages,  ultra-fast  current  exchange  between  two 
collectors  can  be  controlled  by  die  applied  drain  voltage.  Compared  with  CHINT  structure,  the  pro¬ 
posed  RSTLT  employs:  1)  a  recessed  Cl  collector  to  enhance  the  vertical  electric  field  intensity  under 
Cl  for  matched  ou^ut  collector  currents,  2)  unintentionally  doped  channel  and  barrier  layers  for 
reduced  collector-source  and  collector-drain  leakage  currents,  and  3)  a  coUector-up  device  configuration 
to  reduce  parasitic  capacitance  and  to  facUitate  proper  arrangement  of  the  device  terminals. 

In  this  project,  we  employ  a  self-consistent  ensemble  Monte  Carlo  model  to  perform  the  device 
feasibility  study.  Figures  7  and  8  show  the  current-voltage  characteristics  and  transient  current 
response  as  a  function  of  time,  respectively.  A  characteristic  delay  time  of  about  3.0  psec  can  be 
achieved,  which  would  make  the  RSTLT  one  of  the  fastest  switching  devices  reported  to  date.  This 
characteristic  current  switching  (between  two  collector  terminals)  can  be  utilized  to  realize  a  multi-logic 
element  in  a  single  device  structure.  We  have  performed  case  studies  to  demonstrate  that  by  properly 
selecting  device  dimensions  and  the  resistor  load  elements,  collector  terminals  Cl  and  C2  can  function 
as  logic  NOT  and  EQUIVALENT  of  input  terminal  D.  Based  on  these  results,  we  propose  a  RSTLT 
logic  element  (Fig.  9),  which  features  NOR/ AND  functions  (Cl  functions  as  NOR  of  D1  and  D2;  C2 
functions  as  AND  of  D1  and  D2,  respectively)  and  a  comparable  characteristic  delay  time  as  that  of  the 
basic  RSTLT. 

D.  Effects  of  doping  variation  on  electron  transport 

High  performance  electronic  systems  require  devices  with  high  current  drive  capability,  high  tran¬ 
sconductance,  high  breakdown  voltage,  and  low  parasitic  resistance.  Currently,  most  semiconductor 
devices  employ  doping  techniques  that  control  doping  profiles  in  the  vertical  dimension  (e.g.,  uniform 
doping,  modulation  doping,  and  delta-doping),  that  is,  perpendicular  to  the  direction  of  carrier 


Figure  8:  Current  response  for  drain  voltage  switch-off  for  an  AlojCaojAs/GaAs  RSTLT.  The  col¬ 
lector  voltages  are  fixed  at  Vci=Vc2=4.0  V.  The  device  ii^>ut  (Vp)  is  switched  from  4  V  to  1  V  at 
t=1.0  psec. 


Figure  9:  Schematic  cross  section  of  a  RSTLT  NORJAND  logic  element  The  layer  arrangonent  is 
the  same  as  that  in  the  basic  RSTLT  stnicture  shown  in  Fig.  6. 
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transport.  The  rapid  develtqxnent  of  new  technologies  makes  possible  changes  in  doping  profile  in 
both  lateral  and  vertical  (parallel  and  normal  to  carrier  transport,  respectively).  This  added  degree  of 
freedom  for  "doping  engineering"  can  possibly  result  in  novel  and  improved  device  structures.  The 
motivations  of  this  project  are  to  study  the  effects  of  unique  variable  dqnng  profiles  on  electron  tran¬ 
sport  properties  and  to  assess  potential  device  a{^lications. 

We  apply  self-consistent  ensemble  Monte  Carlo  model  to  study  a  GaAs  n'^-n-n'^  test  device. 
Doping  schemes  are  ramp-doping  (rampup  and  ramp-down)  and  spike-doping  (single-spike  and  multi¬ 
spike)  in  the  n  region.  In  order  to  have  meaningful  comparisons  among  these  structures,  we  rq^ly 
conditions  of  identical  minimum  doping  density  and  the  same  total  integrated  dopant  in  the  active  n 
region.  Figure  10  shows  the  current-voltage  characteristics  of  rampup  (from  the  left  n'^-n  junction), 
ramp-down,  and  two  constant  doped  structures.  Figure  1 1  compares  the  current-voltage  characteristics 
for  one-spke  and  three-spike  doped  structures.  It  is  found  for  conditions  of  identical  minimum  doping 
density  and  total  integrated  dopant  in  the  n  region  that  1)  the  most  significant  improvement  can  be 
achieved  by  employing  a  structure  with  a  ramp-down  doping  scheme.  2)  Introducing  n*^  spike(s)  in  the 
active  electrtm  transport  (n)  region  brings  improved  device  performance  compared  to  devices  without 
n"^  spikes.  For  the  one-spike  doped  structure,  electron  transport  is  independent  on  the  location  of  the 
n"^  spike.  However,  less  significant  improvement  is  observed  for  this  structure  due  to  the  presence  of 
retarding  fields  over  a  large  portion  of  the  n  region.  3)  Further  improvement  can  be  achieved  by  split¬ 
ting  one  n'*'  spike  into  several  spikes  in  the  n  region.  The  redistribution  of  electron  density  and  die 
reduced  effect  of  retarding  electric  fields  bring  an  increased  average  electron  velocity  arxi  a  resulting 
higher  current  drive  capability.  These  results  indicate  the  possibility  of  enhanced  device  performances 
for  FET  devices  which  employ  proper  variable  doping  schemes.  Potential  advantages  of  variable  dop¬ 
ing  in  normally-on  FET  device  rqiplications  include  enhanced  current  drive,  reduced  source  resistance, 
and  improved  breakdown  characteristics. 
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Figure  10:  Current-voltage  characteristics  for  ramp-doped  n^-n-n'^  structures.  The  doping  density 
varies  linearly  between  5x10*®  cm"^  and  5x10*^  cm~^  Simulation  results  for  n'*‘-n-n'''  structures  with 
constant  n  doping  densities  of  5x10*®  cm~^  and  1.6x10*^  cm~^  are  included  for  comparison. 
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Figure  11:  Current-voltage  characteristics  for  one-spike,  three-spike,  and  constant  (SxlO'^  cm~^) 
doped  n'^-n-n^  stractures.  The  doping  densities  of  the  n*  spike(s)  and  the  background  are  9.16x10*^ 
cm”^  and  5x10**  cm~^,  respectively.  The  one-spike  structure  has  the  n^-spike  placed  in  the  center  of 
the  n  region. 
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2.2.3  Modeling  of  hot  electrons  in  silicon  devices 

For  the  past  fifteen  years,  the  silicon  MOSFET  industry  has  been  dealing  increasingly  with  prob¬ 
lems  related  to  hot  electron  injection  into  the  oxide  with  production  mode  devices.  For  even  longer, 
laboratories  have  been  measuring  and  modeling  hot  electron  injection  into  the  oxides  with  various  lev¬ 
els  of  success  [25].  Frequently,  the  model  of  choice  has  been  the  "lucky  electron"  model  coupled  with 
a  drift-diffusion  simulator.  While  the  "lucky  electron”  model  has  enjoyed  a  reasonable  amount  of  suc¬ 
cess,  we  have  recently  entered  the  regime  [26]  where  drift-diffusion  simulation  is  no  longer  accurate. 
With  this  reality,  researchers  have  increasingly  been  forced  to  switch  to  hydrodynamic  and  Monte 
Caiio  simulators  to  better  predict  the  behavior  of  devices.  Using  these  models  (especially  Monte 
Cailo),  significantly  more  information  about  the  electron  energy  distribution  is  available,  giving  us  the 
opportunity  to  improve  upon  the  "lucky  electron"  model.  FuU  theoretical  treatment  of  the  entire  hot 
election  injection  condition  poses  many  serious  technical  difficulties  [27]  but  some  first  pass  effions 
[28]  indicate  we  may  have  a  reasonable  hope  for  success.  We  can  then  couple  these  injection  charac¬ 
teristics  with  the  oxide  trapping  characteristics  [29,30]  to  give  us  a  second  generation  model  to  sup¬ 
plant  the  "lucky  electron”  model. 

Within  the  past  year,  we  have  developed  a  state-of-the-art  Monte  Carlo  simulator  of  electrons  for 
silicon  devices.  This  simulator  uses  a  realistic  band  structure  (derived  from  pseudopotential  calcula¬ 
tions)  for  the  best  accuracy  in  simulating  band  effects  on  electron  simulation.  It  also  uses  an  ensemUe 
method  of  modeling  which  allows  the  simulation  to  account  for  both  spatial  and  electron-electron  inter- 
collisional  effects  in  a  self-consistent  maimer.  In  addition,  our  Monte  Carlo  model  includes  all  of 
major  scattering  mechanisms  important  for  silicon  simulation.  These  mechanism:,  are  acoustic  phonon, 
intervalley  phonon,  impaa  ionization,  ionized  impurity  and  electron-electron  scatterings.  This  rather 
accurate  implementation  of  physical  detaib  results  in  a  reasonable  field  vs.  velocity  and  scattering  rates 
which  compare  favorably  with  the  full  numerical  calculations  of  Fischetti  [31]. 
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In  order  to  accurately  model  complex  device  structures,  we  have  also  developed  a  non-uniform 
mesh  Poisson  solver  which  is  re-evaluated  at  frequent  time  intervals  to  reflect  changes  in  the  potential 
distribution  due  to  the  movement  of  the  electrons  in  the  simulation.  This  movement  is  mirrored  in  the 
calculation  of  background  hole  concentrations.  The  carrier  concentrations  remain  self-consistent  by 
updating  the  hole  concentrations  using  the  quasi-Fbimi  level  derived  from  the  Poisson  solution.  At  the 
same  time,  an  efficient  method  has  been  implemented  which  allows  to  transfer  PISCES-modeled  and 
developed  structures  into  one  of  our  Monte  Carlo  runs.  This  approach  will  improve  overall  computing 
time  and  efficiency  significantly  by  providing  a  smart  initial  guess  and.  thus,  will  pennit  us  to  concen¬ 
trate  to  the  regimes  of  sensitivity  in  a  greater  detail. 

Calculations  of  actual  MOSFET  structures  is  now  finishing  the  debugging  stage  and  the  prelim¬ 
inary  data  begin  to  sh  «  .le  importance  and  value  of  our  model  in  the  study  of  hot  electron  effects 
and  the  related  aging  issues  in  these  devices.  Figure  12  demonstrates  that  electrons  experience  a 
significant  velocity  overshoot  region  even  at  reasonably  long  channels  (1  pm)  when  stressed  at  high 
voItag>.s  (Vd=6  V,  Vg=3  V,  V,=Vb=0  V).  Tlie  electron  energy  distribution  tail  in  the  drain  pinch-off 
region  is  also  seen  to  have  a  significant  population  beyond  the  3.1  eV  necessary  to  inject  electrons  into 
the  oxide  as  shown  in  Fig.  13.  Such  details  are  the  essential  information  to  understand  nonlinear  elec¬ 
tron  transport  in  submicron  devices.  We  expea  to  obtain  more  data  in  our  effort  to  develop  a  realistic 
model  for  oxide  trapping  and  the  generation  of  interface  states  at  the  Si/Si02  interface  within  the  next 
few  weeks. 

2.3  Publications  and  Presentations 

During  the  last  year  we  have  made  six  oral  presentations  at  national  and  international  confer¬ 
ences.  Also,  twelve  written  manuscripts  have  been  published  in  the  refereed  literature.  In  addition, 
another  manuscript  has  been  accepted  for  publication  and  six  manuscripts  have  been  sutoiitted  or  are 
currently  in  preparation.  The  following  paragraphs  summarize  the  presentations  and  publications  made 
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under  this  program  during  the  last  year,  and  include  material  not  previously  reported  to  ONR  in  the 
last  annual  report  as  well  as  material  which  was  accepted  for  publication  in  1991  and  has  since  been 
published  in  1992. 
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We  present  results  from  a  theoretical  study  of  the  influence  of  doping  pro^.'c  variations  on 
the  performance  of  delta-doped  AlGaAs/GaAs  high  electron  mobility  transistors 
(HEMTs).  An  ensemble  Monte  Carlo  simulation  coupled  with  a  self-consistent  solution  of 
the  two-dimensional  Poisson  equation  is  used  to  investigate  HEMTs  which  employ 
both  single  and  double  delta-doped  profiles  with  varying  doping  configurations.  The  calculated 
results  reveal  that  single  delta-doped  HEMTs  designed  with  identical  threshold  voltages 
exhibit  improved  device  behavior  when  thinner  delta-doped  layers  with  more  heavily  doped 
concentrations  are  utilized.  For  double  delta-doped  HEMTs  with  an  identical  total 
doping  in  the  AlGaAs  layer,  improved  threshold  voltage  control  is  obtained  as  the  spacing 
between  two  delta-doped  layers  increases.  However,  this  increase  in  spacing  also 
causes  a  degradation  in  transconductance,  cut-off  frequency,  and  switching  time.  As  gate  bias 
increases,  the  dependence  of  device  performance  (or  degradation)  on  the  spacing 
between  doping  planes  becomes  less  pronounced  due  to  the  upward  shift  in  threshold  or 
“onset”  of  parallel  conduction  in  the  AlGaAs  layer. 


I.  INTRODUCTION 

The  high  electron  mobility  transistor  (HEMT)  was 
first  developed  in  1980*  and,  today,  the  AlGaAs/GaAs 
HEMT  remains  as  one  of  the  most  important  heterojunc¬ 
tion  structures  due  to  its  excellent  device  performance,  the 
maturity  of  the  materials  growth  technology,  and  the  ex¬ 
cellent  lattice  match  for  the  AlGaAs/GaAs  material  sys¬ 
tem.  Further  device  improvements  are  needed  in  order  to 
fully  exploit  the  potential  of  AlGaAs/GaAs  HEMT  struc¬ 
tures.  Unfortunately,  improvements  in  “conventional"  uni¬ 
formly  doped  HEMTs  are  limited  by  the  occurrence  of 
persistent  photoconductivity,^  threshold  voltage  shift,^  and 
collapse  of  current-voltage  (/-F)  characteristics^  due 
largely  to  effects  caused  by  DX  centers  and  surface  states. 
The  use  of  low  A1  mole  fraction  can  reduce  these  effects; 
however,  this  t  Jso  reduces  the  conduction-band  edge  dis¬ 
continuity,  resulting  in  decreased  channel  electron  density 
and  degraded  device  performance.  This  inherent  limitation 
can  be  circumvented  by  device  structure  modifications, 
such  as  those  employed  in  delta-doped  (also  referred  to  as 
pulse-doped  or  planar-doped)  AlGaAs/GaAs  HEMT 
structures.  Recently,  delta-doped  HEMTs  have  attracted 
much  attention.  Experiments  show  that  these  devices 
achieve  improved  threshold  control,  high  channel  electron 
density,  r^uced  trapping  effects,  and  high  breakdown 
characteristics.*"'® 

Detailed  investigations  of  delta-doped  HEMTs  are 
necessary  to  reveal  and  understand  the  underlying  device 
physics.  While  extensive  study  of  conventional  uniformly 
doped  HEMTs  has  been  reported,"'"  the  theoretical  in¬ 
vestigation  of  delta-doped  HEMTs  is  rare. The  opera¬ 
tion  of  submicron  delta-doped  HEMTs  involves  nonsta¬ 
tionary  electron  transport  and  real-space  transfer 
effects.""*'  Also,  the  degeneracy  effect,  DX  centers  and 


surface  states  have  significant  influence  on  practical  device 
operation. Therefore,  the  development  of  more  accu¬ 
rate  and  complete  device  models  must  include  these  effects. 
Furthermore,  it  is  important  to  explore  new  approaches 
which  can  lead  to  improvements  in  AlGaAs/GaAs 
HEMTs.  Specifically,  reducing  series  resistance  and  in¬ 
creasing  current  drive  capability  of  the  single  delta-doped 
HEMT  structure  is  of  importance  for  microwave  and 
power  applications.  As  an  alternative,  a  structure  with  a 
second  delta-doped  layer  introduced  in  the  AlGaAs  layer 
(i.e.,  a  double  delta-doped  HEMT)  could  be  used.  This 
structure,  while  unavailable  experimentally,  merits  study 
since  the  introduction  of  a  second  delta-doped  layer  could 
possibly  provide  more  drive  current  and  reduce  source-gate 
and  gate-drain  series  resistance  without  reducing  the  ad¬ 
vantages  provided  by  single  delta  doping.  A  double  delta- 
doped  structure  can  be  realized  by  molecular-beam  epitaxy 
or  metalorganic  chemical  vapor  deposition  technologies. 
However,  a  theoretical  study  can  provide  guidance  for  the 
experimentalists  to  assess  the  performance  of  such  a  device 
and  its  dependence  on  doping  parameters  and  bias  condi¬ 
tions. 

With  these  motivations,  we  employ  physical  device 
models  to  examine  the  effects  of  doping  profile  variations 
on  the  performance  of  two  sets  of  HEMT  devices.  The  first 
set  includes  single  delta-doped  HEMTs  with  varying  dop¬ 
ing  configurations  and  constant  threshold  voltage.  Also, 
for  the  first  time,  an  attempt  is  made  to  analyze  the  per¬ 
formance  of  proposed  double  delta-doped  HEMTs.  lliis 
second  set  of  devices  is  studied  by  varying  the  spacing 
between  two  delta-doped  layers  with  identical  widths  and 
total  doping  concentrations.  Device  performance  is  exam¬ 
ined  in  terms  of  threshold  voltage  control,  transconduc¬ 
tance,  cut-off  frequency,  and  characteristic  switching 
times.  [Since  the  main  purpose  of  this  work  is  to  study  the 
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A  paik-inltfral  Montt  Carlo  method  for  calculation  of  teal  self-energies  of  free 
and  confined  carriers  at  non-teto  temperatures  is  presented.  An  efficient  numer¬ 
ical  method  for  calculating  the  polar-optical  phonon  influtnet  funelionol  hat  been 
developed  and  incorporated  into  a  previously  developed  PIMC  for  calculating  equi- 
librinm  properties  of  uncoupled  carriers.  As  Hlnstrative  examples,  calculations  of 
carrier  real  self-energies  are  performed  for  bulk  GaA$  as  a  function  of  temperature, 
and  for  AlCoAs-CoAr-like  square  quantum  wires  at  77K  and  300K  as  a  function  of 
wire  width.  The  self-energies  obtained  here,  in  part,  represent  new  results  for  these 
systems.  Where  previous  results  are  available,  the  agreement  is  good. 


1.  Introduction 

The  Ftynman  potfi-intcyraf  (FPI)  formulation  of  quan¬ 
tum  mechanics  allows  formal  inclusion  of  carrier  phonon 
coupling  to  all  orders  in  the  coupling  potential  (1,2|.  For 
this  reason  the  FPI  formalism  long  has  been  used  to  ana¬ 
lytically  calculate  carrier  real  self-energies  and  effective 
masses  in  bulk  semiconductors  |3-6|.  However,  exten¬ 
sion  of  the  FPI  formalism  to  analytic  study  of  ultrasmaR 
devices  has  been  hindered  by  approximations  required  to 
model  confined  geometries  |7|.  Previously,  though  im¬ 
proved  results  were  not  obtained,  it  was  demonstrated 
that  carrier  sclf-cnergiea  in  bulk  semiconductors  at  OK 
could  be  calculated  using  numerical  FolA-Inteyrsi  Uontt 
Carlo  (PIMC)  methods  as  weU  (8-10|.  However,  to  the 
best  of  the  authors'  knowledge,  these  calculations  were 
not  extended  to  nonsero  temperatures  or  confined  de- 
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vice  geometries.  Here,  a  PIMC  method  for  calculating 
free  and  confined  carrier  self-energies  at  non-iero  tem¬ 
peratures  is  presented. 

2.  Feynman  Path-Integral  Formalism 

In  the  FPI  formalism,  for  a  coupled  carrier-phonon  sys¬ 
tem,  the  equilibrium  density  matrix  for  the  carrier  av¬ 
eraged  over  the  phonon  coordinates  and  subject  to  a 
Maxwell-Boltsmann  distribution  in  energies  at  temper¬ 
ature  T  is  given  by 

p(f,?iT)  -  j[^xp|-ijf*'‘'|Im*|ri(l)|*.hK|r(l)|ld£j 

(1) 

Here,  the  signifies  an  integral  over  all 

paths  from  P  to  f  that  are  continuous  functions  of  time, 
and  ^rjt)l  it  the  imflutnet  fiinelional  that  accounts  fur 
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Hydrodynamic  electron-transport  model:  Nonparabolic  corrections  to  the  streaming  terms 
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This  paper  presents  a  hydrodynamic  model  suitable  for  studying  hot-electron  transport  in  semicon¬ 
ducting  materials  with  nonparabolic  conduction  bands.  The  model  presented  is  based  upon  a  unique 
derivation  of  the  moments  of  the  Boltzmann-transport  equation  for  the  streaming  (collision-independent) 
terms.  This  derivation  implements  an  efficient  and  compact  mathematical  formalism  appropriate  for 
electrons  under  the  influence  of  high  electric  fields  and  nonstationary  conditions.  The  theoretical  inves¬ 
tigation  also  introduces  a  distributional  form  with  nonparabolic  properties  to  precisely  define  the  result¬ 
ing  nonparabolic  streaming  parameters.  The  final  set  of  model  equations  is  exhibited  in  a  fashion  to 
clearly  show  the  correction  factors  to  the  more  familiar  hydrodynamic  model  applicable  for  the 
constant-effective-mass  case.  In  general,  the  hydrodynamic  (or  conservation)  model  contains  pure  trans¬ 
port  terms  that  are  treated  as  being  independent  of  the  specific  dissipation  mechanisms  and  collision 
terms  to  directly  account  for  the  influence  of  scattering.  Since  the  collision  terms  are  almost  always 
treated  phenomenologically  using  a  relaxation-time  approximation,  our  formulation  of  the  streaming 
terms  should  significantly  improve  the  overall  accuracy  of  the  approach.  In  addition,  this  paper  presents 
the  results  of  an  extensive  investigation  of  the  assumed  ansatz  distribution  and  resulting  nonparabolic- 
model  parameters  using  an  elaborate  Monte  Carlo  model.  The  Monte  Carlo  technique  was  used  to  gen¬ 
erate  comparison  electron  distributions  and  exact  values  for  the  nonparabolic  transport  parameters  for 
stationary  and  nonstationary  electronic  structures.  In  all  cases,  excellent  agreement  was  found  between 
the  Monte  Carlo-calculated  parameters  and  the  derived  nonparabolic-model  terms.  The  Monte  Carlo 
calculations  also  revealed  that  the  ansatz  distribution  used  in  the  derivation  represented  a  significant  im¬ 
provement  over  the  more  familar  displaced  Maxwellian.  Therefore,  this  model  should  prove  very  valu¬ 
able  for  studying  electronic-device  structures  operating  under  high-bias  conditions. 


1.  INTRODUCnON 

During  the  past  two  decades  there  has  been  an  in¬ 
creased  use  of  hydrodynamic  conservation  models  to  in¬ 
vestigate  nonstationary  and  nonequilibrium  electron  dy¬ 
namics  in  submicrometer  semiconductor  devices. 
These  models,  referred  to  by  some”''*  as  the  hydro- 
dynamic  equations  due  to  their  similarity  to  the  Euler 
equations  of  fluid  dynamics'^’'*  used  in  classical  hydro¬ 
dynamics  studies,  are  based  upon  higher  moments  of  the 
Boltzmann-transport  equation  (BTE).  The  popularity  of 
the  hydrodynamic  electron  transport  theory  is  due  to  the 
physical  and  practical  attributes  of  the  a*iproach.  The 
hydrodynamic  models  have  the  capability  to  include  non¬ 
stationary  and  hot-electron  effects'’  and  hence  are  superi¬ 
or  to  simple  drift-diffusion  (extreme  thermal  equilibrium 
approximation^  and  local-electric-field-dependent  case) 
m^els.  While  Monte  Carlo  methods  that  solve  the  full 
BTE  can  easily  incorporate  complicated  band  structures 
and  detailed  scattering  rates,  hydrodynamic  models  re¬ 
quire  much  less  computation  time  to  generate  solutions 
and  possess  macroscopic  terms  which  offer  important 
physical  insight.  Thus  the  hydrodynamic  approach  offers 
much  flexibility  for  future  theoretical  electron  device 
studies. 

Stratton”  introduced  the  general  conservation  or 
momentum-energy  balance  approach  to  investigate  hot- 
electron  transport  in  semiconductors.  The  analysis  per- 
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formed  by  Stratton  utilized  a  spherical  harmonic  expan¬ 
sion  with  the  relaxation-time  approximation  to  define  a 
nonequilibrium  distribution  which  was  only  slightly  an¬ 
isotropic.  This  analysis  replaced  the  nonequilibrium  dis¬ 
tribution  by  a  Maxwellian  distribution  function  in  some 
terms  of  the  BTE.  Thus  some  of  the  streaming  terms  in 
the  final  transport  equations  were  suppressed.  Bio- 
tekjaer^'  extended  the  theory,  using  a  more  general 
analysis  which  retained  all  terms  of  the  moment  equa¬ 
tions,  and  derived  relations  applicable  to  a  two-valley 
semiconductor.  These  previous  investigations  were  appl¬ 
icable  to  semiclassical  transport  and  assumed  homogene¬ 
ous  materials  with  parabolic  conduction  bands. 

Since  the  early  work  of  Stratton  and  of  Blotekjaer,  the 
hydrodynamic  approach  has  been  utilized  extensively  to 
study  various  transport  phenomena  in  many  different 
Si,*’’’”  GaAs,*’**’”  and  heterostructure"’”’"  electron 
devices.  The  particular  material,  device  configuration, 
and  bias  conditions  under  consideration  have  led  to  a 
variety  of  acceptable  simplifications  and  to  a  hierarchy  of 
approximate  electron-transport  models  (for  example,  see 
Sandbom,  Rao,  and  Blakely”  and  references  therein). 
Also,  much  effort  has  been  directed  toward  developing 
numerical  techniques  to  solve  the  hydrodynamic  trans¬ 
port  model  self-consistently  with  Poisson’s  equation  for 
many  semiconductor  structures.  This  has  included 
methods  for  solving  the  problem  with  and  without  time 
dependency  and  in  multiple  space  dimensions”'”  as  well 
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Abstract 

A  new  multi-valley  hydrodynamic  transport  model  has  been  used  to  investigate  electron  transit  time 
through  submieron  GaAs  electron  device  structures.  The  study  was  performed  to  determine  the  feasibility 
of  utilising  the  velocity  overshoot  effect  in  reducing  the  overall  transit  time.  This  new  model,  which 
incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduction  band  valley  (F,  L  and  A') 
and  allows  for  electron  exchange  through  energy  dependent  relaxation  parameters,  is  applied  to  a  novel 
structure  which  reduces  upper  valley  occupancy  and  enhances  device  transit  time.  The  results  of  this 
study  indicate  that  while  velocity  speed-up  is  minimal  in  the  neighborhood  of  the  active  (low-doped) 
region  that  the  velocity  overshoot  can  be  extended  into  the  collector  region  (high-doped)  region  with 
significant  transit  time  enhancement. 


Introduction 

In  recent  years,  a  great  interest  has  developed  in  utilising  velocity  overshoot  in  both  GaAs  homojunction 
and  heterojunction  BJT's  to  achieve  lower  electron  transit  times  and  hence  higher  operating  frequencies)  1 , 
2,  3].  Investigations  of  these  devices  have  shown  the  emitter-collector  transit  time  (Tfc)  and  the  collector 
transit  time  (re)  to  be  dominant  factors  in  determining  the  total  transit  time[4].  This  fact  has  directed 
much  energy  in  the  simulation  and  modeling  communities  toward  designing  and  studying  electronic 
structures  which  would  reduce  these  figures  of  merit  (r,,  and  Te).  In  fact,  there  have  been  many  theoretical 
results  [5,  6,  7|  which  indicate  that  velocity  overshoot  can  be  exploited,  by  select  designs  and  proper  biases, 
to  achieve  very  efficient  high  frequency  devices.  However,  questions  have  been  raised  u  to  whether 
conventional  device  structures  can  exhibit  such  enhanced  effects.  Also,  simulation  evidence  has  been 
reported  to  indicate  that  velocity  overshoot  may  be  over  estimated  in  such  structures[8]. 

In  thu  paper,  a  new  multi- valley  hydrodynamic  transport  model  is  used  to  investigate  the  feasibility  of 
utilising  the  velocity  overshoot  effect  in  reducing  electron  transit  time  through  submicron  GaAs  electronic 
device  structures.  In  previous  work  (9, 10],  a  novel  approach  has  been  presented  to  incorporate  conduction 
band  nonparabolicity  into  the  hydrodynamic  model.  However,  for  this  initial  investigation  with  the  new 
multi-valley  model,  which  incorporates  a  separate  set  of  hydrodynamic  equations  for  each  conduction 
band  (F,  L  and  X)  and  allows  for  electron  exchange  through  energy  dependent  relaxation  parameters, 
the  nonparabolic  correction  terms  will  be  suppressed.  This  was  done  to  simplify  the  analysis  and  to 
provide  a  base  for  future  nonparaboUc  studies.  Additional  features  of  the  model  include  self-consistent 
potentials  and  realistic  boundary  conditions. 
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Abstract.  In  this  paper,  the  results  of  studying  a  submicron  GaAs  ballistic  dirxle, 
using  a  new  multi-valley  (F,  L  and  X)  non-parabolic  hydrodynamic  transport 
model,  are  presented.  Numerical  simulations  indicate  that  accurately  including  the 
effects  of  non-parabolicity  in  the  streaming  terms  and  Ar-space  transfer  in  the 
velocity  and  energy  equations  is  very  important  in  correctly  determining  the 
conductance  of  the  device.  The  existence  and  amount  of  negative  differential 
conductance  was  determined  to  be  strongly  influenced  by  both  of  these  physical 
factors.  Furthermore,  the  sensitivity  of  device  conductance  to  changes  in  the 
thermal  conductivity  is  diminished  significantly  when  non-parabolicity  is 
accurately  incorporated. 


1.  Introduction 

Advanced  macroscopic  transport  models,  based  upon 
moments  of  the  Boltzmann  transport  equation  (bte)  as 
opposed  to  a  full  microscopic  description  (e.g,,  the  Monte 
Carlo  (mc)  approach),  are  valuable  tools  to  efficiently 
investigate  electron  dynamics  in  submicron  semicon¬ 
ductor  structures  [I].  However,  these  models  must  ac¬ 
curately  reflect  non-equilibrium  and  non-stationary 
electron  physics  before  they  can  be  applied  successfully  to 
a  wide  variety  of  device  operating  conditions.  Presently, 
the  majority  of  macroscopic  electron  device  models  treat 
the  multi-valley  problem  indirectly  (effective  single- 
electron  gas  approximation)  and  ignore  the  full  effects  of 
band  non-parabolicity  [2].  When  more  complete  macro¬ 
scopic  models  are  implemented,  they  sometimes  predict 
questionable  results.  For  example,  a  two-valley  (parabo¬ 
lic)  hydrodynamic  model  has  been  used  to  derive 
negative  differential  conductance  in  a  GaAs  fet  while 
corresponding  MC  simulations  indicate  no  such  effect  [3]. 

In  this  paper,  a  new  hydrodynamic  transport  model 
(htm)  suitable  for  studying  hot  electron  transport  in 
multiple  non-parabolic  conduction  bands  is  investigated. 
The  impact  of  (e-space  (intervalley)  transfer  and  band 
non-parabolicity  is  examined  by  studying  variations  of 
this  HTM.  To  determine  the  influence  of  these  effects 
individually,  these  various  forms  of  the  mode)  were 
applied  to  a  submicron  GaAs  ballistic  diode.  Numerical 
simulation  results  indicate  that  including  the  effects  of 
non-parabolicity  in  the  streaming  terms  and  intervalley 
scattering  in  the  velocity  and  energy  equations  strongly 
influences  the  characteristics  of  the  device  in  the  non¬ 
linear  (saturation)  region.  These  results  demonstrate  that 


correctly  including  non-parabolicity  effects  in  the  hydro- 
dynamic  model  reduces  the  amount  of  negative  dif¬ 
ferential  resistance  (nor).  Also,  the  amount  of  nor  was 
found  to  be  less  sensitive  to  changes  in  the  thermal 
conductivity  when  non-parabolicity  is  accurately  in¬ 
corporated.  Therefore,  these  results  verify  the  importance 
of  both  non-parabolicity  and  intervalley  effects  in  the 
HTM. 


2.  The  hydrodynamic  transport  model 

For  these  studies,  a  submicron  GaAs  ballistic  diode  (see 
inset  in  figure  I)  will  be  considered  with  transport  in  each 


Figure  1.  Total  elactron  current  density  against  applied 
electric  potential  for  a  submicron  GaAs  ballistic  diode  (see 
inset).  Results  are  for  variations  in  the  htm  (broken 
curves)  and  a  mc  study  (full  curve).  Here  3  ^  ^  in  the 

HTM. 
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ABSTRACT 


A  novel  real-space  transfer  logic  transistor  (RSTLT)  is 
proposed  based  on  the  real-space  transfer  of  hot  electrons  in 
a  four-terminal  heterojunction  microstmctnre.  Self-consist¬ 
ent  steady-state  and  transient  Monte  Carlo  simulations  dem¬ 
onstrate  that  the  RSTLT  features  both  lope  flexibility  and 
ultra-fast  switching  speed.  Calculated  results  show  that  the 
proposed  RSTLT  realizes  NOT  and  EQUIVALENT  logic 
functions  in  a  single  heterojunction  device  and  a  conser¬ 
vative  estimate  of  the  charaeteiislie  switching  time  is  ~  3 
psec. 


INTRODUCTION 

Real-space  transfer  (RST)  (1]  is  the  thermionic  emis¬ 
sion  of  hot  electrons  (from  one  semiconductor  layer)  over 
heterointerface  barriers  to  different  semiconductor  layers. 
Succeuful  implementations  of  the  RST  concept  have  been 
demonstrated  by  experimental  realizations  of  the  negative 
resistance  field-effect  transistor  (NERFET)  (2|,  charge  in¬ 
jection  transistor  (CHINT)  |3]  in  both  AlGaAs/GaAs  and 
InGaAs/InAlAs  material  systems,  and  real-space  transfer 
transistor  (RSTT)  [4]  uung  strained  InGaAs/AlGaAs/GaAs 
layers.  Promising  performances  have  been  achieved  in  these 
devices,  with  reported  values  of  transconductance  exceed¬ 
ing  2300  mS/mm  (3]  and  unity-current  gain  frequency  of 
60  GHz  |4]  at  room  temperature.  Also,  a  NORAND  func¬ 
tional  element  based  on  a  multi-terminal  CHINT  structure 
has  been  proposed  recently  [5],  demonstrating  the  poten¬ 
tial  logic  applications  of  real-space  transfer  devices.  In  this 
work,  we  propose  a  novel  logic  dremt  element  -  the  red- 
space  transfer  lope  transistor  (RSTLT).  The  operation  of 
this  proposed  microstructnre  is  based  on  the  concept  of  hot 
electron  real-space  transfer  and  the  fact  that  the  spatial 
location  of  electron  RST  is  determined  by  appCed  bias  and 
heteromterfaee  energy  barrier  height.  Self-consistent  ensem¬ 
ble  Monte  Carlo  simulations  are  used  to  theoretically  'build’ 
a  RSTLT  with  flexible  logic  functions  and  ultra-fast  speed. 


RSTLT  STRUCTURE  AND  SIMULATION 
METHOD 


Figure  1  shows  the  proposed  RSTLT  structure.  The  de¬ 
vice  has  two  separate  collectors  (Cl  and  C2)  as  output  ter¬ 
minals  which  are  connected  through  load  elements  (R)  to 
the  positive  power  supply  (V^),  a  drain  (D)  input  termi¬ 
nal  and  a  common  source  (S)  terminal  which  is  at  pound 
potential.  When  D  is  biased  at  logic  lov  (positive  logic), 
channel  electrons  confined  in  the  potential  well  between  the 
channel-barrier  heterointerfaces  will  not  have  enough  energy 
to  overcome  the  energy  barrier  until  reaching  C3  where  they 
are  collected.  In  this  case  Cl  is  logic  high  and  C2  is  logic 
low  (terminal  voltages  of  Cl  and  C2  are  high  and  low,  re¬ 
spectively).  On  the  other  hand,  if  D  is  high,  most  electrons 
gain  suffident  energy  to  overcome  the  barrier  while  passing 
under  Cl.  Therefore,  Cl  becomes  low  and  C2  Kiyh.  Thus, 
the  RSTLT  realizes  the  logic  NOT  (Cl)  and  EQUIVALENT 
(C2)  functions  of  the  input  control  signal  (D),  i.e.,  Cl  s  T5, 
C2  s  D,  as  shown  in  the  inset  of  Fig.  1.  In  order  to  ensure 


figure  1.  Proposed  RSTLT  structure.  The  inset  on  the  top 
right  shows  the  logic  symbol  of  the  structure. 
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Abitraet  •  This  paper  pretenti  an  ealended  Scharfetter-Gummel 
ciirrent-deniitj  diecretiaatlon  formula  eufltcient  for  toWuig  an  ad¬ 
vanced  hydrodynamic  electron  traniport  model.  Thli  Improved  dia- 
cretixation  expreeiion  If  tpecifically  developed  for  a  new  hydrody¬ 
namic  traniport  model  luitable  for  itudylng  hot-electron  traniport 
In  lemiconducting  materiali  with  nonparabolic  conduction  bandi. 
While  thii  extended  formulation  properly  incorporate!  many  ad¬ 
vanced  elTectf  aiiociated  with  hot-electron  traniport,  Iti  primary 
/eatareii  the  accurate  treatment  of  the  eonieetion  term  (or  dtiplacc- 
ment  energy  term)  which  li  preient  in  both  conventional  (l-c.  drifted 
Maxwellian  model)  and  the  advanced  nonparabolie  model  diicnucd 
In  thii  paper.  Our  derivation  overcomei  nonlinear  problemi  intro¬ 
duced  by  thii  convection  term  and  refulti  In  an  Improved  exprenion 
to  approximate  current  deniity. 

I.  INTRODUCTION 

At  preient,  Monte  Carlo  (MC)  methodi  and  hydrodynamic  traniport 
model!  (HTM'i)  are  the  two  moit  popular  approachei  utlliied  to  in- 
veitigate  hot-electron  phenomenon  In  lubmicron  itructurei.  While 
MC  methodi|l|,  which  lolve  the  full  Boltxmann  traniport  erpiation 
(DTE),  can  eaiily  Incorporate  detailed  traniport  phyiici,  thii  ap¬ 
proach  can  become  impractical  in  lomc  initaneei  due  to  the  extcniive 
cempnlationa)  rc<|airemenli.  IITM'i,  which  arc  baaed  upon  macro- 
•copic  eonicrvatlon  equation!  derived  from  momenti  of  the  DTE,  In¬ 
clude  nonitationary  and  hot-eicctron  circcti|2|  while  requiring  much 
left  computation  time  to  generate  lolutioai.  The  phyilcal  and  prac¬ 
tical  attributei  of  IITM'i  have  lead  to  their  extenilve  application.  In 
turn,  advanced  IITM’i  have  been  developed  which  arc  accurate  for 
vorioui  device  itructurei  and  operating  condltioni|3|.  While  IITM’i 
rrpreient  a  coniidcrabic  reduction  in  complexity  (ai  compared  to 
■olving  the  full  DTE),  their  lolutlon  for  arbitrary  device  itructurei 
If  itill  a  very  formidable  motbematlci  problem.  When  IITM’i  ore 
applied  icif-coniiitcntly  with  Poiiioa’i  equation,  one  muit  lolve  a 
let  of  nonlinear  coupled  partial  dliTercnllal  equation!.  Thii  talk  re- 
qiiirei  the  application  of  accurate  and  itabic  numerical  technique!. 
Since  the  lucceii  of  any  numerical  limulotion  li  critically  dependent 
on  convergence  rate  and  numerical  itablilty,  the  choice  of  on  efficient 
diicretiiatlon  icheme  if  very  lmportant|4|. 

It  If  well  known  that  a  central  dllTercncc  approach  ii  not  adequate  for 
diicretlilng  the  current  continuity  equation  due  to  the  development 
of  ofclllationi  or  ‘wipylci”  In  the  lolntlonlel.  To  iniure  numerical 
Itablilty,  Schorfetter  and  Gamincl(6|  Introduced  a  phyilci-baaed  cur¬ 
rent  deniity  formulation  for  the  drlR-dlRaflon  (DO)  model  (l.e.  they 
Illumed  current  deniity  vorlei  weakly  with  poiltlon).  The  lucceif 
the  Seharfctter-Oamniel  (SO)  icheme.  In  OD  model  limulatianl,  hai 
motivated  vorloui  generallaationi  of  thii  diicretliation  technique  for 
the  hydrodynamic  cqaationf|7|.  However,  meat  exteniiuni  of  the 
SO  formula  for  IITM’i  ellher  ignore  convention  termi|4|  or  elae  m- 
lume  e  prforf  that  they  are  defined  In  a  previoai  itep  of  an  Iteration 
procedure(7|.  In  thli  paper,  on  analytical  method  If  pretented  to  de¬ 
velop  an  extended  SG  current  deniity  formula  which  directly  includci 
the  eifecti  of  the  nonlinear  convention  term.  Thli  SG  algorithm  If 
ipecincally  developed  for  uie  In  foivlng  a  new  nonparabofle  HTM. 
Alio,  the  extended  formula  li  compared  to  one  which  Ignoree  energy 
dliplacement  to  liluf  trate  the  Influence  of  thii  nonlinear  effect  on  the 
predicted  current  deniity. 


II.  THE  NONPARABOLIC  MODEL 

Previoiiily,  a  hydrodynamic  model  luitable  for  itudying  boi-eiectron 
traniport  in  lemiconducting  materiali  witii  nonparabolie  conduction 
bandi  wai  preicnted|3|.  Thii  new  model  contain!  advanced  non¬ 
parabolic  itreaming  (colliiion  Independent)  termi  which  were  de¬ 
rived  by  applying  a  unique  let  of  moment  operator!  (l.e.  orderi 
of  nonlinear  k-ipaee  velocity  oi  oppoied  to  momentum)  to  the  col- 
Itfionleii  Boltxmann  traniport  equation.  The  reiulting  lyitem  of 
moment  (balance)  cquotioni  were  cloicd  (terminated)  by  Introduc¬ 
ing  an  empirical  diitribution  function  whirh  accurately  approximate! 
the  propertlei  of  electron!  drifting  in  nonparabolie  conduction  bandi. 
To  treat  Intervally  tranifer  cITecti  that  are  preient  In  multi- valley  lyi- 
temi  (e.g.  GaAi),  a  three-valley  (F,  L  and  X)  eniemble  relaxation 
time  model  hai  been  developed  and  impiemented|8|.  The  full  multi- 
valley  nonparabolie  HTM,  with  a  phenomenological  heat  flow  vector 
term,  coniliti  of  the  traniport  equation! 


^  =  -v..V,w.+P.v,-lv,.(v,.|P,(Km-^qr)  +  (^)^  (3) 


and 

9?  =  (4) 

for  each  conduction  bond  valley  (i  s  F,  L  and  X).  In  the  prcviooi 
model  equationf,  m  b  the  electron  concentration,  V(  h  the  average 
electron  velocity,  w<  if  the  average  electron  energy  and  F  =  yV,d 
b  the  force  on  the  electron.  In  thii  IITM,  the  rciniting  conitita- 
tive  relatiom  arc:  an  average  clfcctlve  mail  ih7(vi,T^)  s  m)(l  -i 
-I-  7.5(a{kprs<)*  -f  aim)|vi|*),  on  cflective  velecity  preffun 
|Fvj|  =  (m)/?h))n<kprM(l|(  on  effective  energy  preiforc  (Pw«|  = 
''ikpTv(l  -f  3aikaTv4  -I-  a{m(|vr|*)(l|,  a  nonparabolie  hoot  low  vec- 
tor  q{  =  koi(m(/iii*)*(kBT«r)*(l  -|-3la((m(/ih))ktT^)n4Xi,  and  the 
nonparabolie  energy  partition  relation  ai,-  =  (3/2)igT,ff(vi,T^)  4 
(*M’’itT»()/2)|v(|*  with  on  clfcetive  temperature  of  T./y(vi,T^)  = 
Ts((I  -f  TatktTmt  4-  a{m)|T(|*).  Theie  previoui  relatioM,  which  ex- 
prcfi  the  corrected  itreaming  termi  aa  functioni  of  v,  w  and  non- 
paroboUc  conduction-band  factor  o  were  verified  for  both  itatiaaary 
and  nonftotiooary  tranaport  in  GoAi  by  two  independent  kleatc 
Carlo  limulation  itndief(3|.  Feotnrea  of  thii  model  which  differ  fieui 
the  conventional  HTM  Include:  (1)  a  noncomtont  average  effectlvt 
man  iIi)(vi,Tb(),  (3)  a  dbtinctiy  dtlferent  vdedly  and  energy  pm- 
lure  (l.e.  IPvil  |i  (Fwrl),  and  (3)  a  flnt-prlnciplei  nonaero  hiat 
flow  vector  q<.  In  thii  model,  qr  repreienti  the  heat  flow  maltlai 
from  the  lymmctiie  eompuncnl  of  the  vclocity.ipace  ehetren  dbtii- 
bution  in  nonparabolie  bandi.  Therefore,  the  WledenMon-Ftaai  phi- 
nomenological  heat  flow  term,  qf  In  Eq.  (4),  b  Introduced  to  accaaal 
for  pofilUe  aa ymmetry  of  the  electron  diitribution.  In  thb  beat  leo 
term,  m  ii  the  near  equllibrlnm  mobfilty  and  Q  b  a  nwdtipllcatlw 
conitant  aaioclated  with  the  particular  material  under  ceneldiritba 
Dbiipatlon  and  Intervalley  tranifer  elfecti  (l.e.  iait  termi  In  Eqa  (1) 

-  (3))  are  Incinded  aalng  an  eneemble  relaxalian  time  model  Idealtal 
In  form  to  the  one  propoaed  by  BlolekJaer|9|.  To  complete  the  madd, 
the  varioui  energy-dependent  raiaxatloa  tlmei  for  a  maltl-vafley  lya 
tom  (GaAf)  have  been  obtained  bom  itationary  MC  catcalatleaiH 
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Abstract  -  The  eifects  of  laterally  varying  doping  profile  on 
device  electron  transport  properties  are  studied  using 
self-consistent  ensemble  Monte  Carlo  simulations.  Significant 
performance  improvement  is  observed  for  N*  -  N  —  N*  struc¬ 
tures  with  ramp-down  doping  and  multi-spike  doping  profiles 
in  the  active  N  region.  The  prospects  for  application  of  later¬ 
ally  varying  doping  schemes  in  field-elTect  transistor  structures 
are  discussed. 


The  utilisation  of  'bandgap  engineering*  has  led  to  a  va¬ 
riety  of  novel  improved  device  structures  in  the  past  decade. 
Successful  integration  of  ‘bandgap  engineering*  with  ‘doping 
engineering*  provides  additional  flexibility  and  new  methods 
for  the  future  development  of  advaneed  semiconductor  devices. 
Currently,  most  semiconductor  devices  employ  techniques  that 
control  doping  profiles  in  the  vertical  dimension  (perpendic¬ 
ular  to  the  direction  of  carrier  transport  using  e.g.,  uniform 
doping,  modulation  doping,  and  delta-doping).  The  rapid  de¬ 
velopment  of  new  technologies  such  at  focused  ion  beams  (FID) 
makes  possible  mask-less  direct  ‘writing*  of  controlled  dopants 
into  semiconductors  |1].  Thus,  the  doping  distribution  can  be 
tailored  in  the  lateral  dimension  (parallel  to  the  direction  of 
carrier  transport).  This  added  degree  of  freedom  for  ‘doping 
engineering*  can  possibly  result  in  novel  device  structures  with 
improved  performance.  The  purpose  of  this  paper  is  to  present 
results  of  a  theoretical  study  of  the  eifects  of  laterally  varying 
doping  profiles  on  device  electron  transport  properties,  which 
supports  this  proposition. 

In  this  work,  self-consistent  ensemble  Monte  Carlo  simula¬ 
tions  are  performed  to  investigate  the  effects  of  laterally  vary¬ 
ing  doping  on  electron  transport.  The  Monte  Carlo  model  em¬ 
ploys  a  three-valley  (P-L-X),  nonparabolic  band  structure  cou¬ 
pled  with  a  solver  for  Poisson *s  equation.  Trajectories  of  a  large 
number  of  sample  electrons  are  traced  in  real  space  and  in  mo¬ 
mentum  space.  Relevant  scattering  mechanisms  (polar  optical, 
equivalent  and  nonequivalent  intcrvalley,  ionised  impurity,  and 
electron-electron)  are  ineluded.  In  addition,  the  effect  of  degen¬ 
eracy  for  P-valley  electrons  is  taken  into  account  by  employing 
an  approximate  Fermi-Dirac  distribution  with  calculated  local 
quasi-Fermi  level  and  electron  temperature  |2|.  Charge  neu¬ 
trality  in  a  portion  adjacent  to  the  electrodes  is  maintained 
throughout  the  simulation,  which  serves  as  the  criterioo  for 
electron  injection.  Material  and  energy  band  parameters  for 
GaAs  are  the  same  as  those  in  Ref.  |3|.  Lattice  temperature 
is  assumed  to  be  300  K.  More  detailed  model  description  can 
be  found  elsewhere  |4|.  This  transport  study  assesses  advan¬ 
tages,  differences,  and  potential  applications  of  different  later¬ 
ally  varying  doping  profiles  on  device  performance.  The  test 
device  is  a  GaAs  Af ♦(0.5/im)  -  Af(0.5/tm)  -  Ar+(0.5pm)  struc¬ 
ture  with  varied  doping  profile  along  the  direction  of  electron 
transport  in  the  N  region.  In  particular,  two  groups  of  laterally 


varying  doping  schemes  are  studied.  They  are:  1)  ramp-doping 
|i.e.,  linearly  increasing  (ramp-up)  or  decreasing  (ramp-down) 
doping  density  within  the  N  region],  and  2)  spike-doping  (i.e., 
introducing  one  or  more  N*  spikes  in  the  N  region).  In  order 
to  have  meaningful  comparisons  among  structures  which  em¬ 
ploy  different  doping  profiles,  we  apply  conditions  of  identical 
minimum  doping  density  and  total  integrated  dopant  in  the  N 
region. 

Any  successful  implementation  of  a  lateral  doping  scheme 
should  result  in  improved  current  drive  capability.  Among  all 
the  devices  studied,  our  simulation  results  reveal  that  the  most 
significant  improvement  can  be  achieved  by  employing  a  struc¬ 
ture  with  a  ramp-down  doping  scheme.  For  the  same  applied 
bias,  the  ramp-down  doped  structure  exhibits  more  than  50  % 
overall  improvement  in  current  drive  capability  compared  to 
that  for  the  structure  with  the  ramp-up  doping  scheme.  This 
is  shown  in  Fig.  1,  where  current- voltage  characteristics  of 
the  structures  with  ramp-up,  ramp-down,  and  constant  mini¬ 
mum  N  doping  profiles  are  given.  The  conduction  band  profile 
(increasing  exponentially  along  the  direction  of  electron  trans¬ 
port)  generated  by  the  ramp-down  doping  scheme  and  the  re¬ 
sultant  accelerating  electric  field  profile  provide  favorable  elec¬ 
tron  transport  conditions  in  which  electrons  accelerate  more 
efliciently.  High  average  electron  velocity,  low  average  electron 
energy,  and  high  F-valley  electron  occupancy  are  obtained  for 
the  ramp-down  doped  structure.  These  results  indicate  that, 
on  average,  more  electrons  in  the  ramp-down  doped  structure 
accelerate  while  they  stay  in  the  F  valley.  Results  from  spike- 
doped  structures  show  that  introducing  an  N*  spike  in  the 
active  electron  transport  (N)  region  brings  improved  device 


Figure  I.  Current- voltage  characteristics  for  the  ramp-doped 
N*  —  N  —  N*  structures.  Doping  density  varies  liaearly  be¬ 
tween  S  X  10'*em**  and  5  x  I0'’em"*.  Result  for  Af*  -  JV  — 
N*  structure  with  constant  (minimum)  N  doping  density  of 
5  X  I0'*em~*  is  included  for  comparison. 
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AlInAs-GalnAs  hemts 

U.K.  Mishra*,  L.M.  Jelloian+,  M.  Lui'*',  M.  Thompson+/  S.  E.  Rosenbaum+ 
and  K.W.  Kim# 

•University  of  California,  Santa  Barbara,  CA.  93106 
+Hughes  Research  Laboratories,  Malibu,  CA.  90265 
#N.C.  State  University,  Raleigh,  N.C.  27695 

1.  Introduction: 

Enhancement  in  the  high  frequency  performance  of  transistors  is  dependent 
both  on  improved  materials  properties  and  smaller  transit  distances.  This  has 
resulted  in  improved  performance  as  devices  have  progressed  from  MESFETs 
(Wang  etal)  to  GaAs  based  pseudomorphic  HEMTs(Tan  etal)  to  InP  based 
lattice-matchedfMishra  etal,  Chao  etal)  and  pseudomorphic(Thompson  etal) 
HEMTs.  The  GalnAs  based  devices  have  demonstrated  the  highest  speeds 
because  of  a  combination  of  high  electron  mobility,  peak  velocity  and  sheet 
charge  density  in  the  GalnAs  channel  of  an  AlInAs-GalnAs  HEMT.  The 
reduction  in  gate  length  to  50  nm  has  increased  the  extrinsic  fr  of  the  HEMT 
to  292  GHz  at  room  temperaturefThompson  etal)  but  the  fmax  of  the  HEMT 
has  been  restricted  by  the  rapid  increase  in  the  output  conductance  at  the 
small  gate  length.  The  output  conductance  is  caused  by  a  combination  of 
substrate  injection  and  channel  length  modulation.  Attempts  to  reduce  the 
substrate  injection  by  using  a  combination  of  p-doped  and  wide  band  gap 
buffer  layers  have  had  limited  success.  In  this  study,  we  investigate  the  effects 
of  doping  the  channel  of  an  AlInAs-GalnAs  modulation  doped  transistor  n 
and  p  type  and  evaluate  the  effect  on  the  DC  and  RF  characteristics  of  devices 
with  0.25  pm  gate  length.  Devices  with  1  pm  gate  length  were  also  studied  to 
determine  the  gate  length  dependence.  Hie  motivation  is  to  evaluate  the 
effect  of  the  two  doping  types  on  the  distribution  of  electric  field  in  the 
channel  and  study  its  effect  on  output  conductance  and  electron  transport. 


GalnAs  Conlact  Layw  7  nm 


Undopad  AHnAs  Layer  20  nm 

Layer 

.  (6ElB/cm3:12  nm) 


GalnAs  Channel  Undoped  or  doped  9 


n(1E17rcm3) 
p(1e17  or  3E17/cm3) 


AllnAs  Buffer  250  nm 


InP  Substrate 


Figure  1:  Schematic  of  the  Layer  Structures  Used  in  the  Study. 
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Abstract— Key  design  parameters  for  delta-doped  GaAs 
MESFET’s,  such  as  delta-doping  profile,  top  layer  background 
doping  density,  and  scaling  of  lateral  feature  size,  are  investi¬ 
gated  using  a  two-dimensional  numerical  simulation.  A  three- 
region  (delta-doped  conducting  channel,  top  layer,  and  sub¬ 
strate)  velocity-field  relation  is  implemented  in  the  model  as 
appropriate  for  the  particular  device  structure  which  is  simu¬ 
lated.  Simulation  results  show  excellent  agreement  with  our 
fabricated  0.5-pm  gate-length  delta-doped  GaAs  MESFET’s 
based  on  atomic  lay er  epitaxy  material .  An  extrinsic  transconduc- 
tance  of  370  mS/mm  and  a  drain-source  current  of  270 
mA/mm  are  obtained  for  typical  devices,  and  the  maximum 
transconductance  is  as  high  as  400  mS/mm,  which  are  the  best 
dc  results  yet  reported  for  0.5-pm  gate-length  delta-doped  GaAs 
MESFET’s.  Considerations  of  design  and  optimization  are  dis¬ 
cussed  in  terms  of  threshold  voltage  sensitivity,  transconduc¬ 
tance,  current  drive  capability,  and  cutoff  frequency,  based  on 
both  simulation  and  experimental  results. 


I.  Introduction 

The  drive  to  achieve  high-speed  and  high-fre¬ 
quency  electronic  systems  has  led  to  continued  efforts 
to  develop  new  device  structures.  One  veiy  promising 
structure  is  a  field-effect  transistor  (FET)  which  incorpo¬ 
rates  the  delta-doping  technique.  Advanced  delta-doped 
FET  structures  such  as  GaAs  MESFET’s  [l}-[4], 
HEMT’s  on  GaAs  [5]-[7]  and  InP  substrates  [8],  as  well 
as  pseudomorphic  HEMT’s  [9]  have  been  demonstrated. 
The  advantages  of  delta  doping  in  these  device  structures 
include  high  current  drive  capability,  high  transconduc¬ 
tance,  improved  threshold  voltage  control,  and  improved 
breakdown  characteristics. 

Delta-doped  GaAs  MESFET’s  based  on  materials 
grown  by  MBE  and  MOCVD  materials  have  been  inves¬ 
tigated  experimentally  by  several  research  groups  [l]-[4]. 
While  the  results  are  promising,  they  also  show  a  wide 
diversity  and  fall  short  of  theoretical  predictions  [1].  The 
reported  maximum  transconductance  of  O.S-pm  gate- 
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length  delta-doped  MESFET’s,  for  example,  ranges  from 
75  mS/mm  [1]  to  140  mS/nun  [2],  to  the  400  mS/mm 
transconductance  reported  in  this  paper.  Schubert  et  al. 
[1]  studied  properties  of  delta  doping  in  GaAs  and  pre¬ 
dicted  that  the  maximum  intrinsic  transconductance  of  a 
delta-doped  MESFET  could  be  higher  than  that  of  a  com¬ 
parable  AlGaAs/GaAs  HEMT  in  the  gate  length  range  of 
0.1  to  2.0  nm.  Differences  between  experimental  results 
and  theoretical  predictions  may  stem  from  unknown  ma¬ 
terial  properties  and  differences  or  uncertainty  in  process¬ 
ing  conditions.  On  the  other  hand,  design  and  optimiza¬ 
tion  of  such  devices  play  an  important  role  in  guiding  the 
improvement  of  device  performance.  It  is  of  practical  im¬ 
portance  to  identify  the  key  parameters  and  the  tradeoffs 
needed  for  optimized  device  performance  in  specific  ap¬ 
plications. 

The  intent  of  this  paper  is  to  present  a  generalized  anal¬ 
ysis  of  delta-doped  submicrometer  MESFET’s.  We  use  a 
two-dimensional  drift-difiusion  model  to  simulate  the  in¬ 
fluence  of  key  design  parameters  on  delta-doped  MES¬ 
FET  performance.  Experimental  0.5-jim  gate-length 
delta-doped  MESFET’s  grown  by  atomic  layer  epitaxy 
(ALE)  have  also  been  fabricated.  Excellent  agreement  be¬ 
tween  simulation  and  experimental  measurement  has  been 
obtained.  An  extrinsic  transconductance  of  370  mS/mm 
and  a  drain-source  current  of  270  mA/mm  are  measured 
for  typical  devices,  while  the  maximum  transconductance 
is  as  high  as  400  mS/mm.  These  are  the  best  experimen¬ 
tal  results  yet  reported  for  0.S-;xm  gate-length  delta-doped 
GaAs  MESFET’s.  Based  on  the  agreement  between  sim¬ 
ulation  results  and  experimental  results  for  O.S-fim  de¬ 
vices,  we  then  analyze  performance  dependence  on  key 
device  parameters  and  outline  some  design  considerations 
and  tradeoffs  for  achieving  improved  delta-doped  MES¬ 
FET  operation. 

n.  Device  Fabrication 

A  schematic  illustration  of  a  delta-doped  GaAs  MES¬ 
FET  is  shown  in  Fig.  1(a).  A  very  narrow  doping  profile 
(Fig.  1(b))  is  sandwiched  between  two  undoped  (or  rela¬ 
tively  low-doped)  layers.  The  device  strucmre  and  the 
doping  profile  in  the  figures  are  representative  of  those 
employ^  in  the  experimental  device,  although  the  actual 
device  fabricated  for  measurement  in  this  paper  contains 
a  recessed  gate.  The  doping  profile  used  in  the  simulation 
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Ensemble  Monte  Carlo  simulations  are  employed 
in  order  to  explore  the  feasibility  of  a  novel  real-space  transfer 
logic  transistor  (RSTLT)  structure.  The  operational  principles 
of  the  proposed  RSTLT  are  based  on  the  concept  of  hot  electron 
real-space  transfer  (RST),  including  the  fact  that  the  spatial 
location  of  electron  RST  is  determined  by  applied  bias  and  het¬ 
erointerface  energy  barrier  height  in  a  multiterminal  hetero¬ 
junction  microstructure.  The  results  of  two-dimensional,  self- 
consistent  steady-state  and  transient  simulations  demonstrate 
that  the  proposed  RSTLT  features  ultrafast  current  switching 
which  can  be  used  to  realize  NOT/EQUIVALENT  logic  func¬ 
tions  in  a  single  heterojunction  device.  The  logic  operation  is 
easily  extended  to  NOR/AND  functions.  A  conservative  esti¬ 
mate  of  the  characteristic  delay  time  for  current  switching  is 
—  3  ps  in  the  proposed  RSTLT  structure. 


I.  Introduction 

HERMIONIC  emission  of  hot  electrons  (from  one 
semiconductor  layer)  over  heterointerface  barrierfs) 
to  different  semiconductor  layers  in  the  presence  of  high 
electric  fields  is  termed  real-space  transfer  (RST)  [1]. 
Successful  implementations  of  the  RST  concept  have  been 
demonstrated  by  experimental  realizations  of  the  nega¬ 
tive  resistance  field-effect  transistor  (NERFET)  [2]-[4], 
the  charge  injection  transistor  (CHINT)  in  both  Al- 
GaAs/GaAs  [5],  [6]  and  InGaAs/InAlAs  [7],  [8]  mate¬ 
rial  systems,  and  the  real-space  transfer  transistor  (RSTT) 
using  strained  InGaAs/AlGaAs/GaAs  [9]  heterostiuc- 
tures.  Promising  performances  in  both  dc  and  microwave 
applications  have  been  achieved  in  these  devices,  with  re- 
poned  transconductance  values  exceeding  2300  mS /mm 
[6]  and  unity-current  gain  frequency  of  60  GHz  [9]  at 
room  temperature.  Recently,  a  CHINT  logic  element  to 
realize  nor/and  functions  [10]  has  been  proposed,  dem¬ 
onstrating  the  potential  logic  applications  of  real-space 
transfer  devices.  Theoretical  investigations  of  electron 
transport  properties  and  related  device  physics  in  real- 
space  transfer  structures  have  also  been  performed  [Ill- 
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[18].  These  studies  provide  a  foundation  for  device  struc¬ 
tures  with  novel  operating  principles. 

It  is  worth  noting  that  in  a  real-space  transfer  device 
(such  as  the  CHINT),  hot  electron  RST  is  a  localized  pro¬ 
cess  in  which  the  spatial  occurrence  of  RST  is  determined 
by  factors  such  as  the  device  configuration,  heterojunc- 
tion  material  parameters,  and  operating  bias  conditions. 
This  phenomenon  has  been  discussed  by  interpreting  ex¬ 
perimental  I-V  characteristics  of  CHINT  structures  using 
positive  feedback  initiated  by  a  local  ‘hot  spot’  in  the  de¬ 
vice  channel  [8].  Also,  Monte  Carlo  analysis  of  three- 
terminal  RST  devices  has  indicated  the  localized  nature 
of  electron  RST  [13].  Recently,  we  have  proposed  a  novel 
real-space  transfer  logic  transistor  (RSTLT)  [19]  based  on 
ensemble  Monte  Carlo  studies  of  hot  electron  transport  in 
a  four-terminal  heterojunction  microstructure.  The  oper¬ 
ation  of  the  proposed  RSTLT  is  based  on  the  fact  that  the 
spatial  location  of  electron  RST  is  determined  primarily 
by  applied  bias  and  heterointerface  energy  barrier  height. 
The  dc  and  transient  current  switching  characteristics  of 
the  proposed  RSTLT  have  been  discussed  briefly  [19] 
(here  we  define  current  switching  as  the  current  exchange 
among  device  terminals  while  transient  switching  times 
for  the  intrinsic  device  are  denoted  as  characteristic  delay 
times).  We  have  demonstrated  the  the  potential  logic  ap¬ 
plications  based  on  multiple  terminal  current  switching 
characteristics  of  the  RSTLT. 

In  this  work,  detailed  theoretical  investigations  using 
self-cor  'stent,  two-dimensional  ensemble  Monte  Carlo 
simulatii  s  are  presented  which  focus  on  exploring  the 
feasibility  of  the  RSTLT  structure  in  digital  applications. 
The  proposed  structure  differs  from  the  CHINT  logic  ele¬ 
ment  [10],  i.e.,  the  RSTLT  employs  a  different  device 
configuration  and  terminal  arrangement  to  achieve  shorter 
transient  delay  time  while  maintaining  logic  flexibility. 
This  provides  another  attractive  approach  to  realize  ultra¬ 
fast  multilogic  functions  in  a  single  heterojunction  de¬ 
vice.  This  computer  experimental  study  is  used  to  theo¬ 
retically  “build”  a  device  which  allows  the  logic 
flexibility  and  fast  speed  of  available  real-space  transfer 
devices  to  be  explored.  Electrical  characteristics  and  re¬ 
lated  device  physics  of  the  proposed  RSTLT  structures 
are  examined  in  terms  of  material  parameters,  device  fea¬ 
ture  size,  and  bias  conditions.  In  the  next  section,  we  de- 
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Theoretical  results  of  electron  transport  in  n-type  lno52Alo4gAs  are  presented.  The  transport 
properties  of  this  important  semiconductor  were  obtained  using  the  Monte  Carlo  method.  In 
particular,  velocity-electric  field  characteristics  for  different  temperatures  and  doping 
concentrations  in  bulk  In^  ]2Alo.4gAs  are  calculated  for  the  first  time.  Physical  parameters  for 
100  32^10.48^5  (which  is  lattice-matched  to  InP  and  Ino33Gao.47As)  were  obtained  based  on 
interpolation  of  available  experimental  and  theoretical  results  for  InAs,  AlAs,  and 
Ino  7sA1o  2)As.  Our  study  suggests  that  Ino.32Alo.4gAs  has  electron  transport  properties  which  are 
comparable  to  and  complimentary  with  those  of  other  materials  lattice-matched  to  InP. 


Heterojunction  devices  based  on  the  lattice-matched 
Ino.32Alo.4gAs/Ino33Gao47As/lnP  material  system  have 
been  the  focus  of  extensive  experimental  studies  due  to 
their  promising  potential  for  digital,  microwave,  and  opti¬ 
cal  applications.  Impressive  device  performance  has  been 
achieved  for  the  high  electron  mobility  transistor 
(HEMT),'  charge  injection  transistor  (CHINT),^  and 
semiconductor  laser^  fabricated  from  this  material  system. 
However,  certain  material-related  issues  remain  which 
merit  careful  study  if  potential  applications  of  these  devices 
are  to  be  fully  exploited.  As  an  example,  high  output  con¬ 
ductance  has  been  repeatedly  observed  in  deep-submicron 
( ^0.3  /um)  channel-length  HEMTs  using  the 
I%5:Alo.4gAs/Ino,33Gao.47As/InP  material  system. This 
performance  degradation  has  been  suggested  to  be  the  re¬ 
sult  of:  ( 1 )  electron  conduction  in  the  lno.32Alo4gAs  buffer 
layer;  (2)  deep  level  traps  in  the  top  Ino.32Alo4gAs  donor 
supply  layer;  and  (3)  weak  impact  ionization  in  the 
lno.33Gao47As  channel  layer.  Studies  on  the  CHINT  have 
also  shown  that  electron  transit  time  in  the  In^  32AIo4gAs 
collector  barrier  is  the  main  factor  that  determines  the 
total  device  intrinsic  delay  time.^  Th^  observations  indi¬ 
cate  the  importance  of  achieving  a  better  understanding  of 
electron  transport  in  Ino.32Alo4gAs.  Unfortunately,  there  is 
very  little  available  information  on  material  parameters  or 
electron  transport  properties  for  InQ32Alo4gAs.  This  is  a 
noticeable  void  since  there  have  been  extensive  experimen¬ 
tal  and  theoretical  studies  of  InP  and  Ing  33Gao47As,  and 
data  on  Ino32Alo4gAs  are  required  in  order  to  perform 
thorough  and  more  rigorous  device  simulations  in  the 
Ino.S2Alo48As/Ino  33Gao.47As/InP  material  system. 

In  this  letter,  we  employ  a  Monte  Carlo  model  to  study 
electron  velocity-electric  field  characteristics  in  n-type 
Ino32Alo4gAs  for  different  temperatures  and  doping  con¬ 
centrations.  The  set  of  physical  parameters  for 
Ing  32Alo.4gAs  required  for  this  study  has  been  obtained 
based  on  linear  interpolation  of  available  data  for  the  InAs 
and  AlAs  binary  components,  while  using  the  few  available 
(theoretical)  material  properties  of  Ino73Alo2}As  as  a 
guide.  The  material  parameters  which  we  provide  for 
Ino  32Alo.4gAs  <‘"<1  results  of  the  fundamental  electron 
transport  study  will  serve  to  promote  further  experimenta¬ 
tion  and  characterization  of  Ioq  32Aio.4gAs-related  devices. 


The  Monte  Carlo  method  used  in  this  work  comprises 
a  model  for  the  motion  of  sample  electrons  in  a  three  valley 
( r -L-X)  semiconductor  with  nonparabolic  analytical  band 
structure  under  the  influence  of  a  uniform  applied  electric 
field.^’^  Relevant  scattering  mechanisms  are  polar  optical 
phonon  scattering,  equivalent  and  nonequivalent  interval¬ 
ley  scattering,  ionized  impurity  scattering,  alloy  scattering, 
and  impact  ionization.  In  the  present  work,  a  single 
longitudinal-optical  (LO)  phonon  mode  is  used.  There  is  a 
possibility  that  two  LO  phonon  modes  due  to  the  InAs-like 
and  AlAs-like  vibrations  could  be  active  in  Ino32Alo4gAs, 
since  it  is  known  that  this  condition  does  exist  in  some 
ternary  materials.  However,  this  is  a  subtle  point,  and  de¬ 
tailed  properties  of  the  phonon  spectra  in  this  material  are 
not  sufficiently  known  at  present  to  justify  the  use  of  two 
LO  phonons  in  the  model.  The  use  of  a  single  LO  phonon 
with  an  energy  of  39.5  meV,  which  was  interpolated  from 
local  mode  energies  of  29.6  meV  for  InAs  and  50.1  meV  for 
AlAs,  should  not  detract  from  the  main  conclusions  of  this 
work. 

The  ionized  impurity  scattering  rate  used  in  our  sim¬ 
ulations  is  based  on  Ridley’s  statistical  scattering  model.^ 
In  this  model,  the  scattering  rate  in  terms  of  electron  wave 
vector  k  is  given  as 


Ti„p(k)  d  «;,(k)TBH(k)/j’ 


(1) 


where  i^gfk)  is  the  electron  group  velocity  and 
</[=(27rA^dop)  ~'^^]  '5  Ihe  average  distance  between  impu¬ 
rities.  The  electron  density  is  assumed  to  be  equal  to  the 
doping  concentration  (A^op).  The  scattering  rate  l/rgH. 
which  is  derived  from  the  Brooks-Herring  model,^  was 
calculated  by  following  Ruch  and  Fawcett. 

The  alloy  scattering  rate  is  obtained  using^ 


1 


’*’illoy 


(x(l-x)]y(e) 


1/2 


dy{€) 

de 


n|At/|^5(a), 
(2) 


where  AC/  is  the  alloy  scattering  potential,  x  is  the  alloy 
composition,  11  is  the  primitive  cell  volume,  and  y(e) 
=e(  1  -l-ae)  for  the  nonparabolic  conduction  band.  We  use 
the  electron  affinities  of  AlAs  and  InAs  to  predict  a  value 
for  AC/.*  Here,  S(a)  is  an  energy-dependent  parameter 
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